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Abstract 

Warped scenarios offer an appealing solution to the hierarchy problem. We con- 
sider a non-trivial deformation of the basic Randall-Sundrum framework that has a 
KK-parity symmetry. This leads to a stable particle beyond the Standard Model, that 
is generically expected to be the first KK-parity odd excitation of the radion field. We 
consider the viability of the KK-radion as a DM candidate in the context of thermal 
and non-thermal production in the early universe. In the thermal case, the KK-radion 
can account for the observed DM density when the radion decay constant is in the 
natural multi-TeV range. We also explore the effects of coannihilations with the first 
KK excitation of the RH top, as well as the effects of radion-Higgs mixing, which im- 
ply mixing between the KK-radion and a KK-Higgs (both being KK-parity odd). The 
non-thermal scenario, with a high radion decay constant, can also lead to a viable sce- 
nario provided the reheat temperature and the radion decay constant take appropriate 
values, although the reheat temperature should not be much higher than the TeV scale. 
Direct detection is found to be feasible if the DM has a small (KK-parity odd) Higgs 
admixture. Indirect detection via a photon signal from the galactic center is an inter- 
esting possibility, while the positron and neutrino fluxes from KK-radion annihilations 
are expected to be rather small. Colliders can probe characteristic aspects of the DM 
sector of warped scenarios with KK-parity, such as the degeneracy between the radion 
and the KK-radion (DM) modes. 



1 Introduction 



Models of electroweak symmetry breaking (EWSB) that contain a stable particle closely tied 
to the EWSB sector open a window of complementarity between astrophysical observations, 
DM searches, and collider experiments. Not surprisingly, such a prospect has attracted 
considerable attention in the context of supersymmetry, little Higgs models with T-parity [1] 
or models with Universal Extra Dimensions (UED) [2J. 

The latter case is interesting in that the reason the DM candidate is stable is closely tied 
to the geometry/compactification of the extra dimensional space [HIS]- In the 5D case, this 
is simply due to a Z 2 reflection symmetry (usually called KK-parity) about the center of 
the compactified fifth dimension. The viability of the lightest KK-parity odd particle (the 
LKP) as a dark matter (DM) candidate was first studied in [51 [6]. For a review, see [7]. An 
important assumption in the previous construction is that the extra dimensional spacetime 
is flat. On the other hand, since the seminal work of [E], the presence of a non-negligible 
spacetime curvature (although with vanishing 4D effective cosmological constant), has been 
understood to have rather important and interesting consequences. Unfortunately, in typical 
implementations, no stable new particles appear in such warped scenarios, although it is 
simple to build models with a stable DM candidate by simply imposing a Z 2 symmetry 
(see e.g. [H]). There have also been proposals where the DM and EWSB sectors share non- 
trivial properties [TOl ITT1 IT2] , within the framework of warped compactifications, as well as 
studies in "warped throats" that appear in type IIB string theory on Calabi-Yau manifolds 
where a suitable DM candidate can be obtained [13]. Also, Ref. p3j has considered the 
possible Sommerfeld enhancements arising from radion exchange for GUT-motivated DM 
candidates in warped scenarios, while stable candidates arising from a dark sector in warped 
compactifications have been proposed in Ref. [T5] . 

It is nevertheless possible, and not difficult, to build warped scenarios with a KK-parity 
symmetry, similar to UED models. The first attempt was done in [16j, where it was proposed 
to glue together two copies of the RS model along either the UV or IR branes. This reference 
considered the possibility that the first KK-parity odd excitation of the Z-gauge boson could 
be the DM candidate. Nevertheless, it is well known that the lightest state beyond the SM 
in such constructions is expected to be the radion field [17J. If this field has a KK-parity 
odd excitation, it is more likely that this "KK-radion" will be the LKP. A concrete model 
that stabilizes the interbrane distance, and hence gives mass to the radion field, while at the 
same time leading to the existence of a KK-parity symmetry, was proposed in Ref. [18J (see 
[19] for an alternative proposal). Indeed, it is found that generically the LKP is the first 
KK-radion mode. 

In this work, we explore the viability of the KK-radion as a DM candidate. Given that 
the radion is such an intrinsic feature of extra-dimensional scenarios, having it constitute 
most of the matter density of the universe would be a rather pleasing possibility. The KK- 
radion mass is found to be parametrically lighter than the scale of the typical Kaluza-Klein 
(KK) resonances, which is set by k^ = k e Re~ A<yL \ where k e g is the curvature at the (IR) 
boundaries of the metric background ds 2 = e~ 2A ^r) fJ , u dx tJ 'dx 1 ' — dy 2 . In models where the 
KK scale is of order TeV, so that they can be tested at the LHC, the LKP mass is expected 
to be around the EW scale. Like the radion field itself, it couples to the SM particles (and 
to other KK excitations) via higher-dimension operators, suppressed by the "radion decay 
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constant" [IB] 



Ar = J X ~k eS , (1) 

where M 5 is the 5D (reduced) Planck mass. We note that, from a theoretical point of view, 
A r can take on a large range of values by varying the ratio M 5 /k e g, while keeping k eS ~ TeV. 
In fact, it was pointed out in Ref. [TS] that varying the ratio M 5 /k eS allows to interpolate 
between RS-like "strong warping" scenarios (when k e ^ ~ M5) and "UED-like" scenarios 
(when k e g <C M5). Although the former case may be considered more natural, we will also 
consider the second possibility, since it gives an interesting deformation of UED models, that 
shares properties of pure UED and the RS model. Furthermore, this is also an interesting 
limit from the point of view of the KK-radion relic density. 

As is well-known, obtaining the correct relic abundance of an EW scale thermal relic 
requires that its interactions be of typical EW size. Thus, unless A r is sufficiently close to 
the EW scale, the KK-radion interactions would be too weak, the KK-radion would freeze-out 
rather early in the history of the universe, and the relic density following from the thermal 
freeze-out paradigm would be too large and lead to the "overclosure" of the universe. As 
we will show, however, the scenario is nevertheless viable over a wide range of A r . First, it 
turns out that the self-annihilation cross section for KK-radions can be sufficiently large even 
when A r is in the multi-TeV range. Second, the presence of other relatively light colored 
resonances can lead to a sufficiently strong depletion of KK-radions through coannihilation 
effects, and make the scenario viable for A r as large as O(10 A TeV). Mixing with a KK 
excitation of the Higgs can have a similar effect. For much larger values of the radion decay 
constant, it is possible to imagine that the KK-radions never reached thermal equilibrium 
in the early universe. This is the case in the "UED-like" (or "small- warping" ) scenario 
mentioned above. It is then necessary that the non-thermal production of KK-radions be 
sufficiently suppressed. This restricts the range of allowed A r for a given reheat temperature, 
Tr. For instance, if Tr ~ 1 TeV, one must have A r ~ 10 15 GeV. 

The outline of this paper is as follows. In Section [2j we summarize the relevant properties 
of the scenario proposed in Ref. [18J. In Section [3j we comment on the expected finite 
temperature properties of our setup. In Section |4j we consider the KK-radion relic density in 
the thermal freeze-out scenario, while we consider the effects of coannihilations in Section [5] 
The case of non-thermal production and the related constraints are treated in Section [6] 
We also comment on the prospects for direct and indirect DM searches for KK-radion DM, 
as well as on the collider signals that can help in distinguishing this scenario from other 
possibilities. We conclude in Section [8j We provide two appendices where we discuss the 
KK-radion couplings and the effects of EWSB (Appendix [A]), and where we summarize the 
relevant cross sections and decay rates (Appendix [B]). 

2 Models with Warped KK-Parity 

We start by summarizing the salient properties of the scenario presented in Ref. [IB]. This 
is a 5D model with the fifth dimension compactified to an interval, and parameterized by 
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y G [—L,L\. The size of the extra dimension is stabilized at tree-level by a bulk (GW) 
scalar field, $, similar to the Goldberger-Wise proposal [17J. However, it is assumed that 
the gravitational backreaction of the stabilizing scalar is non-negligible, and that it is the 
main source for a non-trivial warping of the 5D spacetime. When the scalar potential 
satisfies V{— <&) = V($), and allows for a kink-like background solution ($(?/)) = <f){y), 
with 4>(—y) = —<fi(y), the associated warp factor, e~ A( ~ y \ is found to be even about y — 0, 
i.e. one finds A(— y) = A(y). The non-trivial kink configuration is assumed to be chosen 
over the trivial vacuum expectation value (VEV), <f>(y) = 0, by appropriate potentials on the 
infrared (IR) boundaries, at y — ±L. The background scalar energy density is concentrated 
around y = 0, thus leading to a dynamically generated "UV brane" at the origin. In general, 
this is a "fat" brane with thickness governed by the mass parameter appearing in the bulk 
scalar potential, although there is a well defined limit where the UV brane is thin and the 
background has a constant, negative curvature. In this limit, one recovers the IR-UV-IR 
configuration considered in Ref. [To] . 

For the purpose of the present work, the most important consequence of the above setup 
is the existence of a discrete Z-i symmetry (KK-parity), under which the KK modes whose 
wavef unctions are even about y = are assigned KK-parity +1, while those whose wavef unc- 
tions are odd about y = are assigned KK-parity —1. The resulting KK-parity symmetry 
is a good quantum number, and the lightest KK-parity odd excitation (the LKP) is exactly 
stable. 

It is also important to note that a generic feature of the spectrum is the relative degener- 
acy between the KK-parity even and KK-parity odd states of a given 5D field. The reason is 
that what differentiates even from odd states are the properties of the corresponding wave- 
functions, f n (y), at the origin: KK-parity even states satisfy /4(0) = 0, while KK-parity odd 
states satisfy / n (0) = 0. But when these wavefunctions are localized near the IR bound- 
aries, so that their values are (sometimes exponentially) suppressed at the origin, only a 
small deformation is required to convert an even wavefunction into an odd one. Therefore, 
the associated KK masses are very close. This observation is particularly interesting when 
applied to IR localized (would-be) zero-modes. We highlight three important examples. 

First, the lightest physical scalar excitation of the 5D metric/bulk scalar system is a 
KK-parity even field that is typically referred to as the "radion" . Its wavefunction is given 
by F (y) « e 2[A{y) - A{L)] . This state has a mass that is parametrically lighter than the typical 
KK scale. From the argument above, the first KK-parity odd scalar excitation has a mass 
that is exponentially close to the radion mass (at tree-level). This "KK-radion" is expected 
to be the LKP, and shares properties closely connected to those of the radion field. We will 
refer to it either as r' of r_. 

Second, since a solution to the hierarchy problem requires the Higgs field to be highly 
localized near the IR boundaries, one finds light KK-parity even and KK-parity odd SU (2) 
Higgs doublets. The KK-parity even doublet is fully responsible for EWSB (giving rise 
to a SM-like Higgs and the would-be Goldstone modes eaten by the W ± and Z). The 
KK-parity odd Higgs is "inert" in the sense that it does not participate in EWSB, but it 
can be parametrically lighter than the KK scale since, were it not for EWSB, it would be 
exponentially degenerate with the KK-parity even Higgs doublet, at tree-level. However, we 
will assume that it is somewhat heavier than the KK-radion, since it is a generic scalar whose 
mass is not protected by any symmetry. Annihilation of KK-radions into the KK-parity even 
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Higgs degrees of freedom will be found to play a relevant role in setting the DM relic density. 

Third, fermion zero-modes can be localized either near the "UV brane" (at the center 
of the fifth dimension) or near the IR boundaries. An appealing possibility is to use this 
localization properties to understand the observed fermion mass hierarchies [201 !2"Tj . Thus, 
assuming that the 5D Yukawa couplings are all of the same order, the lighter SM fermions 
(which include the RH tau and bottom) are localized in the UV region, while the top quark is 
localized closer to the IR branes. In particular, one might expect that the RH top is localized 
close to the IR boundaries, and therefore its first KK-parity odd excitation is expected to 
be the lightest of the fermionic resonances. How light this state is depends on how closely 
localized to the IR branes the RH top is. We will refer to this KK-parity odd particle as t' . 
The next-to- lightest KK-parity odd particle (NLKP) is either t' or the inert Higgs doublet 
above. 

It should be pointed out that in the anarchic scenarios, where the 5D Yukawa couplings 
are sufficiently large to generate the top mass, with the (t, b) SU (2)-doublet as close to 
the first two generations as possible (to minimize potentially large corrections to the Zb^L 
coupling), the effects from EWSB are non-negligible. We take these into account by diag- 
onalizing the KK mass matrix in the presence of the EWSB contributions. A consequence 
of this is that the lightest KK-parity odd fermion resonances are pushed downward in mass 
("level repulsion" in the absence of a zero- mode for the KK-parity odd subtower). As a 
result, we find that the t' cannot be heavier than about a factor of 4 above the radion 
mass, and is typically much closer to m r ~ m r /. However, one might also wish to consider 
non-anarchic scenarios where the Yukawa couplings themselves are hierarchical [22, 23j, and 
EWSB effects are much less important. 

To summarize, the states lighter than the KK scale, k e g, are the radion, its first KK-parity 
odd mode (which is stable), a light KK-parity odd t' not much heavier than the previous 
two, and an additional Higgs doublet, on top of the SM field content. Our main objective 
in this work will be to study the viability of this scenario from the point of view of DM 
constraints. If the KK radion is indeed the LKP, these are strongly dependent on the radion 
decay constant given in Eq. ([I]), and we will explore several possibilities. 

3 High- Temperature Phase Transitions 

Before presenting the analysis of the KK-radion relic density, we comment on the expected 
high-temperature properties of the present scenario. In the well-studied case of the RS back- 
ground with the radion stabilized via a Goldberger-Wise scalar, the transition from a high- 
temperature "deconfined phase" to a low-temperature "confined" one has been discussed 
in Refs. [21 [25j EQ EH E3 EH ED] • The terminology arises from the AdS/CFT correspon- 
dence [3TI [32] , where the phase transition maps into a deconfinement /confinement phase 
transition in the dual, large N gauge theory (this phase transition is strongly first-order). 
The analysis has been performed on the gravity side, and in the limit that the backreaction 
of the stabilizing GW scalar is small, by comparing two finite-temperature solutions obey- 
ing the same boundary conditions on the UV brane (however, Ref. [29] considered the case 
where the backreaction is large). These are: the Euclidean RS background, with an IR brane 
and Euclidean time (as measured on the UV brane) compactified on [0,/3 = 1/T]; and the 
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AdS-Schwarzschild (AdS-S) solution, where the IR brane is replaced by a black-hole horizon. 
The free energy associated with the RS solution vanishes as a result of the requirement that 
the effective 4D cosmological constant be zero at the minimum of the radion potential. The 
free energy of the AdS-S solution is then given by F dec ~ E Q — cN 2 T 4 , where E ~ m^A 2 . is 
a vacuum energy contribution, and the second term corresponds to a plasma of deconfmed 
"gluons" with N 2 degrees of freedom at temperature T (c is an order one number). Recalling 
that under the correspondence, A r ~ Nk, it follows that the transition temperature is of 
order T c ~ (m 2 . A; 2 ) 1 / 4 , where k is the warped-down curvature scale. 

Our setup differs from the above RS scenario in two ways: there is no "UV boundary" , 
and the backreaction of the stabilizing dynamics is not a perturbation (but rather the major 
source of the warping). Also, the analog of the AdS-S solution is not known in this setup, 
although it is plausible that such a solution -with black-hole horizons symmetrically located 
with respect to y = 0- exists, and has similar properties to the one mentioned above (recall 
that there is a limit where our scenario reduces to two copies of RS glued together at the UV 
brane). By analogy, one might expect that the scaling relations of the previous paragraph 
also hold in the present case, with k — > k e g. In addition, below T c , the rate of true-vacuum 
bubble nucleation is expected to be suppressed by e~ 0( - N For moderate values of N, i.e. A r 
near the KK scale, one can then expect that, as the temperature drops below T c , there will 
be a period of supercooling, and an associated inflationary era while the universe remains 
trapped in the false vacuum, as discussed in (2U |27] for the RS case. Eventually, bubbles 
nucleate and the prompt decay of the radion reheats the universe to a temperature of order 
TeV. From then on the history of the universe is standard, and the KK-radion relic density 
is set by the dynamics of states that are parametrically lighter than the KK scale, and are 
in thermal equilibrium. The result is unambiguous and independent of the precise reheat 
temperature, T R . 

However, we will also be interested in scenarios with small warping (fc e ff "C M5), hence 
A r ^> k C R. Since this corresponds to the very large N limit, we expect that the rate of bubble 
nucleation will be extremely suppressed. To avoid an empty universe, we therefore need to 
assume that the highest reheat temperature never exceeded T c , and that the universe was 
in the "confined" phase throughout its history. Note that in our scenario m r ~ k eS /^k cS L 
is at the EW scale for any A r {^] Therefore, T c is expected to be at the TeV scale, and is 
independent of N. Nevertheless, to account for possible uncertainties in the above analysis 
that may allow T c to be slightly larger, we will also study, at least formally, the case in which 
the reheat temperature is as large as an order of magnitude above k e g (roughly the 4D cutoff 
of the KK theory). With this in mind, we proceed to describe several scenarios that can lead 
to KK-radion DM. 

4 The KK-Radion as a WIMP 

We first consider the r' relic density in the standard thermal freeze-out scenario. As men- 
tioned in the introduction, all the interactions of the KK-radion, like those of the radion itself, 
are non-renormalizable and controlled by the radion decay constant, A r . An important ex- 

1 In terms of the dual 4D picture language, the explicit breaking of the conformal symmetry is in no sense 
small at the scale where this symmetry is spontaneously broken. 
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Figure 1: Self-annihilation of r' into fermion pairs. The u-channel diagram is understood. 



ception arises from a possible mixing with the Higgs, and we will consider this case separately 
at the end of this section. It is therefore a question whether these non-renormalizable in- 
teractions are sufficiently strong to deplete the thermal r' density down to acceptable levels, 
which will set an upper bound on A r . A possible loophole arises when A r is so large that 
it is possible to imagine that the KK-radions never reached thermal equilibrium in the first 
place. We postpone the discussion of such a scenario to Section [6j Assuming then that the 
r' are in thermal equilibrium at a high temperature, the self-annihilation into SM particles 
are responsible for maintaining a thermal number distribution when the temperature falls 
below the KK-radion mass, m r >. 

The KK-radions can annihilate into fermion pairs, as shown in Fig. [T] The annihilation 
into massless gauge bosons are suppressed at tree-level by (m r '/myj) 4 , where m VJ is the 



mass of a gauge KK mode, as well as by a small (volume suppressed) coupling (see Eq. (68) 



of Appendix B.l in Ref. [IE]). We shall therefore ignore these channels in the following. The 
annihilation into the W and Z gauge bosons proceed through their longitudinal polarizations, 
and can be roughly taken into account through the annihilation into Higgses shown in Fig. [2] 
(this approximation will be sufficient for our purpose). 

We give the full annihilation cross sections in the non-relativistic limit in Appendix [B] 
[Eqs. (63)-(67)]. Note that, as indicated in Fig. [TJ the annihilation into fermion pairs involves 
a sum over the fermion KK tower in the t and li-channels. Including, for illustration, a single 
state f\ the total annihilation cross section into fermion pairs can be well approximated by 



N c m 2 r , 
4vrA4 



t 1 " 4") 



3/2 



1 + e 



~-f,r> 



+ 6 /,r< 



,(2) 



where N c = 3 (N c = 1) for quarks (leptons), €fy = mf/m r >, e-jy = mj/m r >, and G 
characterizes the r'-fjf vertex, while Gnoo characterizes the r'r'-ff vertex, n It would seem 



from the second term in the square brackets, that the cross section increases like mj for 
heavy KK states. However, in the same limit, the coupling constants Gijq decrease even 
faster, and one can check that the sum over j is typically dominated by the first couple of 
terms. The effects from EWSB, which are sizable for the top tower in anarchic scenarios, are 
important in enhancing the above effective couplings. In particular, we find that for the ti 
channel, Gnoo ~ ~6 due to numerical factors characteristic of the two-KK-radion couplings 
[e.g. the factor of 4 in A^ jfc , given after Eq. (58)], as well as due to EWSB mixing with the 



2 To obtain this simplified expression we took g^k <§; (assuming tz, is nearly flat), and we also neglected 



gijQ (for the top, this one is not small, but ends up giving a subdominant contribution 

QLR 



We further took 

Gyo (= Gi 3 o) and defined Gnoo = Giwo = ^uoo (^ ne notation is defined in Appendix [b|) . 
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Figure 2: Self-annihilation of r' into Higgs pairs. The u-channel diagram is understood. 



KK-tower. For the first top KK resonance, although Gno ~ 0.3, the enhancement due to the 
larger KK mass makes the second term in the square brackets of Eq. (|2| non-negligible (but 
still subdominant). It is also important to note that GijQ is proportional to EWSB. In fact, 
the full expression ([2| is proportional to mj, although this fact is somewhat hidden inside 
Gijo for the second term in square brackets. Therefore, the annihilation rate of KK-radions 
into the light families is negligible compared to the annihilation into top quarks. 

The annihilation of KK-radions into Higgses proceeds through the diagrams in Fig. [2] 
In the limit of strong IR localization of the Higgs field, there is a KK-parity even Higgs 
doublet, H + , that is responsible for EWSB, and a KK-parity odd doublet, H_, that does 
not participate in EWSB, but that can be parametrically lighter than the KK scale (see 
comments in Section [2]). We assume in this work that this KK-parity odd Higgs doublet is 
heavier than the KK-radion (at the end of this section we comment on the possibility that 
the CP-even component of the "inert" Higgs doublet is the LKP). In principle, there can 
be heavy Higgs resonances contributing in the t and it-channels, as for the annihilation into 
fermions, but these resonances decouple in the limit that the Higgs is strongly IR localized. 
For simplicity, we will restrict ourselves to this limit, which will be sufficient to illustrate our 
point. If the SM-like Higgs is lighter than the KK-radion, while the inert Higgs doublet is 



heavier, we have [see Eq. (67) in Appendix IB 



N H m 2 r > fo\ 

Va r ' r '^HH ~ 7-7— ■ [d) 



r 



where, for large m r >, N H = 4 takes into account the SM-like Higgs as well as the would-be 
Goldstone modes, and therefore also annihilation into W , Z. 

In order to estimate the total KK-radion annihilation cross section, we assume that 
the RH top tower is localized near the IR boundaries, that the LH top is approximately 
uniform along the extra dimension, and that the lighter families are localized near y = 
(the UV region). For concreteness, and using the parameterization discussed in Ref. [32] to 
characterize fermion localization, which is analogous to the one used in the pure Randall- 
Sundrum AdSs background, we take c tR = —0.2 and c tL = 0.52 (these reproduce the top 
quark mass for k e g ~ 1.2 TeV, after including EWSB effects). As explained before, we can 
neglect annihilations into the light families, but need to take into account the annihilation 
into the Higgs and massive gauge bosons. For reference, in the "strong warping" benchmark 
scenario defined in Ref.[T8] and used in this section, the branching fractions are: BR(rV — > 
tt) ~ 0.70 and BR(rV — > HH) 0.30, where H denotes the four real d.o.f. in the Higgs 
doublet. It is expected that these are distributed as BR(rV — > hh) BR(rV — > ZZ) « 
|BR(rV ->• W+W~) ~ 7.5%. 
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Figure 3: Contours of constant Q r /h 2 in the m r /-A r plane, in models of flavor anarchy, 
where EWSB effects are significant. We use the "strong warping" benchmark relation m t i w 
2.377V, so that coannihilations give a small contribution, and take a SM-like Higgs mass 
rrih ~ 150 GeV, while the KK-parity odd doublet has mass tuh = 3m r r. 

Given the total thermally averaged annihilation cross section, (voyv), the r' relic density 
is given by [33J 

„ , 2 1.04 x 10 9 GeV" 1 x F 

Vr'h 2 « — r , (4) 

Mpi yfg*(V0- rlrl ) 

where is the effective number of relativistic degrees of freedom at the time of r' decoupling, 
Mpi = 1.22 x 10 19 GeV is the Planck mass, and xp = m r //Tp gives the freeze-out temperature 
by solving for 



(5) 

X=Xp 

where one can take ( ~ 0.5, and in this case g = 1. 

In the left panel of Fig. [3j we show contours of constant Q r /h 2 in the m r >-A r plane, 
including the previously described annihilation channels. Since we are assuming a flavor 
anarchy scenario, we take into account EWSB effects, which are non-negligible in that case. 
We take a SM-like Higgs mass of rrih = 150 GeV, while the mass of the KK-parity odd Higgs 
doublet is taken as m# = 3m r r. As mentioned before, we take into account the annihilation 
into gauge bosons by using a multiplicative factor iV# = 4 for the KK-parity even Higgs 
channel. We vary m r > by scaling k e g, i.e. keeping all mass ratios in the spectrum fixed. 
However, we force the SM fermions masses to remain fixed at their observed values as k e g is 
changed. The dimensionless couplings arising from overlap integrals are fairly insensitive to 
k eS , and are kept fixed to the ones corresponding to the "strong warping scenario" defined 
in Ref. [18]. The radion decay constant can be controlled via the ratio M 5 /k e s, which takes 



x F = log|c(C + 2)v / 45^( 



g \ m r 'Mp(va r / r /) 



2vr 



/g*x 



values between about 1.0 and 2.4 in the plot. We note that in this example m t i/m T i 2.3 
with 77V ~ 0.22 k e s, so that the NLKP is sufficiently heavier than r' that coannihilation 
effects are negligible. As can be seen from the figure, the correct relic density of Q r rh 2 rs 0.1 
can be obtained for natural values of the parameters. In the "strong warping" benchmark 
scenario with k e g = 1.2 TeV, we have m r i = 270 GeV and A r = 4.4 TeV, which leads to 
Q r >h 2 ~ 0.24. Reducing A r to 3.6 TeV gives agreement with the WMAP constraint. One 
finds xp ~ 22, so that the freeze-out temperature is about 10 GeV. 

Thus, we see that in spite of the non-renormalizable interactions of the radion, the 
"WIMP miracle" can be operative. This is due to sizable radion couplings to fermion pairs 
(associated with EWSB) and can be qualitatively understood from our expression (§. Tak- 
ing fees = 1-2 TeV, m r i = 0.25 k e s, rnj = m t i ps 0.5& e ff, Gnoo ~ —6 and Guo ~ 0.3, one 
estimates (oyvv/c) ~ 0.84 pb x (3 TeV/A r ) 4 , where we recall that WMAP requires a cross 
section of about 0.8 pb. 

One should also keep in mind that the radion interactions can be enhanced in the presence 
of the operator RH^H, where R is the Ricci scalar. Such a term leads to kinetic mixing 
between the radion and CP-even, KK-parity even Higgs [34, 35J, as well as between the KK 
radion and the CP-even, KK-parity odd Higgs. As a result, renormalizable interactions to 
the SM particles are induced (although the mixing angles are formally suppressed by v/A r ). 
The presence of such a mixing offers the best chance at direct detection of KK-radion DM, 
and we defer the associated technical discussion to Section [7} However, here we comment 
on the effects for the relic density computation in the presence of such mixing. The two 
lightest KK-parity odd mass eigenstates are in general admixtures of r' = r_ and h = /t_, 
which we call <\>-,l and 0-,//. The subscripts L (H) refer to the lighter (heavier) of the two 
mass eigenstates. The DM candidate is 0dm — 4>-,l- Only when the mixing is small, or the 
unperturbed m r i is much smaller than , can this state be thought as mostly KK-radion 
(our assumption in the previous analysis). There is a similar mixing in the (r, h) = (r + , h + ) 
sector, thus defining mass eigenstates +i l and Both of these KK-parity even states 

contribute in the final state when they are lighter than the DM particle, although it may 
happen that only one is sufficiently light to contribute. The Feynman diagrams of Figs. [T] 
and [2] are generalized to those in Fig. |4j Note the presence of an s-channel diagram, which 
can lead to a resonant enhancement if any one of the KK-parity even masses is close to 
about twice moM = J71-,l- The final states include also the Goldstone modes, although 
those processes are not changed by the mixing with the radion tower. 

The mixing with the Higgs field can significantly enhance the DM annihilation cross 
section, and therefore allow for larger values of A r . As a point of reference, in the limit that 
irih- < m r' and A r — > oo, the DM becomes pure h-, and its annihilation cross section, which 
is dominated by the tt channel, reduces to 

N c x t R tj ml 

where X tRt \ is the coupling of h- to the RH top and the first KK resonance of the LH 
top. Fl This coupling is relatively large, and therefore the result is dominated by the t- (and 
it—) channel diagrams of the annihilation into fermions. With X tnt i ~ 2 (which includes 



3 Thc third top resonance has a larger coupling that compensates for its larger mass, so that it ends up 



9 



E 

3 </>DM 



'DM S 



+ 




!>DM \ 



+ 



/ <i>+,L 
' H 



"DM 



H 



H 



4>+.L 
H 



'DM S 



H 



+ x 

/ V 

i ' s j 

0DM ' S <P+,L 



0DM \ ✓ 0+,L 

S V H 

/ ^+- L V 

' H \ 

0DM 7 X 4>+,L 

H 



Figure 4: Self-annihilation of 0dm — in the general case where the radion/Higgs and 
KK-radion/ inert Higgs mix. The subscripts L (H) refer to the lighter (heavier) state, while 
the subscripts + and — refer to KK-parity. These generalize the diagrams of Figs. [I] and\^ 



EWSB effects), and taking m t > ~ rn DM w (but still assuming that m t ' > rn^), we find 
that the observed relic density is obtained for moM ~ 30 TeV. Such large masses are not 
well motivated in the present scenario, and it follows that for my l _ ~ 1 TeV, a too large 
h- component will deplete very effectively the final DM density. We conclude that the DM 
should be mostly KK-radion, although a small h_ component that is interesting from the 
point of view of direct DM searches is allowed. 



5 Coannihilations with t' 

We have seen in the previous subsection that the r' self-annihilation cross section can be 
large enough to avoid the overproduction of KK-radions, provided the radion decay constant 
is sufficiently low. Does this mean that models with much larger radion decay constants are 
ruled out? Here we emphasize the point made in Section [2] that the first KK top excitation 
is expected to be relatively close to m r >. The reason is two-fold: first, the SU(2) doublet 
and singlet top fields cannot be too far from the IR brane, or else m t cannot be reproduced. 
This means that the first KK-parity odd excitation of either t^ or tn must be parametrically 
lighter than the KK scale. Second, when EWSB effects are included this state is pushed 
down even further due to the level repulsion of the second KK-parity odd state. 

If it turns out that m t i is relatively degenerate with the KK-radion (but heavier), coan- 
nihilations with the strongly interacting t' can become dominant. In fact, when A r is large, 

giving a comparable contribution to the t' exchange. The heavier KK modes give a much smaller effect. 
The Higgs annihilation channels involve the Higgs quartic coupling, A. For a heavy Higgs (i.e. h+), these 
channels, if kinematically open, can give a contribution comparable to the tt one. 

4 Assuming that the 5D top Yukawa coupling saturates the NDA [35] estimate of Y 5 d/L < 3/(47r), and 
taking both top zero-modes as far away as possible from the IR boundaries (ct R ~ 0.45 and Ct L ~ 0.48), 
we find that before EWSB the lightest KK-parity odd eigenvalue is about ~ 2.31£; ff, which goes down to 
~ 0.89 k e g when EWSB is taken into account. In the "strong warping" scenario considered in Ref. [18] . with 
Ct R ~ —0.2 and Ct L ~ 0.52, we have me « 0.5fc e ff- The KK-radion mass is m r i ~ 0.22fc c ff. 
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Figure 5: Self-annihilation oft' into gluon and fermion pairs via the QCD interactions. 



the final r' relic density can be completely controlled by the t' self-annihilation cross sec- 
tion, and hence by QCD (see Fig. [5]). In this case, the coannihilation channel shown 
in the left diagram of Fig. [6] gives a negligible contribution to the effective annihilation 
cross section. However, it is important that the crossed diagram, r'g — > i't, although 
suppressed by (mt/A r ) 2 , can be effective in depleting the KK-radion number density, and 
thus maintain the r' in equilibrium until the time the t' decouple. Indeed, for my ~ m r i 
and at the temperature associated with the decoupling of t', i.e. m t '/T F ~ 25, one has 
n g (vo(r'g i't)) /n t ,(va(t'i' gg,qq)) ~ {T F /m t f/ 2 e m t' /TF x (m t /A r ) 2 . Thus, even with 
A r as large as ~ 10 4 TeV, the r' number density maintains an equilibrium distribution by 
processes such as r'g — > i't (or the crossed reactions t'g — > r't) until t' freeze-out. 

The final relic density, after the KK tops have decayed into the KK-radions, is determined 
by an effective cross section 



Ocff 



0~ r / r /^.ff t HH + 9t' e XA (1 + A) 3/2 (a r i t >^ gt + O r i V ^ g t + 0- r 't'^ht + Vr't>-*hi) 



vV 



9cs 



+ gle~ 2xA (l + A) 3 ( Vt'P^qq + OTft^gg 

1 + 9f e~ xA (l + A) 3//2 . 



(7) 



where g t > = 2 x iV c is the number of degrees of freedom associated with t' (and t'). Here 
x = m r //T, while 



A 



ruf — m r i 



parameterizes the degree of degeneracy between r' and t' . Although it is straightforward 
to include them, we neglect EW processes involving annihilations through a Z or a photon, 
since they give a small contribution compared to those involving gluons. 



We give the expressions for a r > t i^ gt and a r i t >-^ht in Eqs. (69) and (70) of Appendix |B| 
respectively. Here we focus on the case where A r is sufficiently large that all processes 
suppressed by powers of A r can be ignored (except, implicitly, for the purpose of maintaining 
thermal equilibrium, as discussed above). As a result, and provided that t! and r' are 
sufficiently close in mass, the final relic density is largely determined by the QCD-controlled 
cross sections 



vovf/ 



-v,/<; 



2na 2 
9m 2 
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27m 2 ' 



(9) 



11 



r' 1 — * — t 



r' 1 — > — t 

g r wvr L ^ — t' 



Figure 6: The left diagram (as well as a diagram where the gluon is replaced by a Higgs) can 
play a role in the co- annihilation scenario. The right diagram is an example of the processes 
responsible for maintaining the r' in thermal equilibrium until t' freeze-out. 

where, for simplicity, we neglected the SM quark masses. It is straightforward to keep the full 
dependence on m t , and we will do so in the numerical analysis, but the effects are relatively 



minor (much smaller than the errors involved in the approximation of Eq. (11) below to the 
full integration of the Boltzmann equations). Similarly, the velocity-dependent terms in all 
the cross sections make a negligible difference. In summary, we have the thermally averaged 
effective cross section 

1727m 2 (1 + A) 3 

{VCTcS) ~ 3m 2 [e* A + 6(l + A)3/ 2 ] 2 ' 1 ' 

The relic density is then approximately given by [37J 

„ l2 1.04 x 10 9 GeV" 1 x F 

Vr'h 2 « — - , 11 

Mpi yfglla 

where = 86.25 is the effective number of relativistic degrees of freedom at the time of t' 
decoupling, the annihilations after freeze-out are taken into account by 

la « XpTdx^-, (12) 

J Xp X 

and xf = m t '/Tp gives the freeze-out temperature by solving Eq. ^ with the replacements 
(voyv) -> (vcr cS ) and g -> g cS . 

Imposing the 7-year WMAP constraint, Qh 2 = 0.1109 ± 0.0112 [38J, determines A as a 
function of m t i (or m r i). In Fig. [7| we show this constraint, as a function of m r i. We see that 
the relic density can be accommodated for a wide range of r' masses, provided the required 
degeneracy with t' is present. We also show in the same figure the freeze-out temperature, 
as parameterized by xp. 



6 The KK-Radion as a Non-Thermal Relic 

We now wish to consider a qualitatively different scenario: one where the KK-radion never 
reaches thermal equilibrium. This can happen when the radion decay constant is very large. 
As was explained in Section [2j this can be realized in scenarios that share properties of the 
Randall-Sundrum models and UED's. Hence they can be motivated outside the context of 
the DM relic density. We assume that the mixing with the Higgs is negligible. 
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Figure 7: WMAP constraint on the degree of degeneracy, A, as defined in Eq. (FM), as a 
function of m r i (red band, right axis). We also show the associated xp (blue band, left axis). 
We assume that A r = 100 TeV ; which is large enough for the relic density to be controlled 
by the QCD t' annihilation cross section, as well as A. 



As discussed in Section |3j when A r ^> k c g we must assume that the reheat temperature 
is lower than the critical temperature for the deconfinement /confinement phase transition. 
Otherwise, the universe would effectively be trapped eternally in the false vacuum. Even 
if a bubble of true vacuum nucleated, the ensuing universe would be dominated by radion 
oscillations until their decay, which would happen rather late. The resulting reheat temper- 
ature may not be high enough to allow for a successful BBN. If instead, after a period of 
inflation, the inflaton decay reheats the universe only to Tr < T c we can assume that the 
universe was always in the "confined" phase. Due to the superweak couplings of the radion 
field, we can assume that it remains at its zero-temperature minimum, so that no energy is 
stored in the radion tower. Note also that, since the radion mass is of order the EW scale, it 
is easy to suppress the deSitter fluctuations of the radion field during the inflationary epoch. 
Thus, radion oscillations are expected to play no role, unlike in the case of UEDs with a 
very light radion Provided the reheat temperature is not much lower than the EW 

scale, KK-radions can be produced in the scatterings or decays of KK-parity odd particles 
in the plasma, but at a slow rate due to the suppressed couplings. Nevertheless, it is well 
known that such production can produce a significant number of superweakly interacting 
particles [40] . which can translate into a bound on the reheat temperature Tr if overclosure 
of the universe is to be avoided. 

As we will see, this requires the reheat temperature to be relatively low, around the 
cutoff of the extra-dimensional theory or lower, even if we overlook the upper bound set 
by T c . In addition, it bounds A r from below to ensure that the KK-radions are sufficiently 
weakly coupled. An upper bound on A r can also be obtained by requiring that the late 
decays of the NLKPs or the non-thermally produced radions do not conflict with Big-Bang 
nucleosynthesis (BBN) constraints. 
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6.1 Non-Thermal Production in Scatterings and Decays 

Our point of departure is the Boltzmann equation 

d ^ + 3Hn r , = C, (13) 

where H is the Hubble constant and n r > is the number density of r'. The collision operator on 
the r.h.s. contains several r' production channels that are our main concern in this section. 
However, since we are assuming that the r' are produced very inefficiently, so that they never 
reach thermal equilibrium, we neglect the inverse processes that would tend to deplete the 
r' number density. This assumption will be checked self-consistently at the end. 

The production channels can be divided into scattering and decay processes. The former 
include scatterings of thermal KK fermions off the plasma, such as gP — > r't k , where P and 
t h are KK excitations of the top quark (including EWSB effects, these are admixtures of the 
doublet and singlet SU(2)l top towers). The towers of the other SM fermions (the "light 
families" ) can also contribute, and will be included in our analysis. In addition, at sufficiently 
high temperatures, the KK excitations of the gluons, g n , can also be active and lead to r' 
production, e.g. via g n gi — > r'g k or g n q^ — > r'q k . Nevertheless, as we will see below, the 
scattering processes represent a subdominant source of thermal r' production, and therefore 
we will be content with obtaining an estimate for the r' production via scattering based on 
the processes with a zero-mode gluon and KK fermions given above. For the same reason we 
also do not consider scattering processes involving the weak gauge bosons]^] We should point 
out that the production of higher KK resonances of the radion can also eventually produce 
r'. However, as was emphasized in [TS], the interactions of such states are controlled by a 
larger decay constant [A" ~ x n A r , where m n = x n k e g e~ A ^ is the n-th KK-radion mass], 
and therefore this contribution is negligible. 

KK-radions can also be produced in (rare) decays of KK tops, P — > r't k , as well as of 
the KK excitations of the light families, and from gauge boson decays, V-[ — > r'V k , with 
Vfj, = G^,W^, Z^, A^. We also include the decays of the inert Higgs doublet, h_ — > h + r', 
a — > G°r' and H ± — > G ± r\ where h + is the SM-like Higgs and G°, G are the would-be 
Nambu-Goldstone bosons eaten by the Z and W ± . Again, it is safe to neglect the production 
of heavier KK-radion resonances. 

During a radiation dominated epoch the Boltzmann equation can be written in terms of 
x = m t i /T as 

where h r '(x) = x 3 n r >(x) and H(x) = H(m t /)x~ 2 = 1.66 gl ' 2 m 2 /M P i x x~ 2 , with g* the 
effective number of relativistic degrees of freedom at temperature T. Note that we choose 
to measure the "time variable" x in units of the t' mass, m? (we are assuming that t' is 



the NLKP). Eq. (14) can be integrated immediately from xr = mt< /Tr, where Tr is the 
reheating temperature, up to any given x, thus obtaining n r >(x) and thereby n r '(x). 

Note that most of the heavy particles involved remain in thermal equilibrium until they 
disappear (by decaying) from the bath: the KK fermion resonances can decay via P — > t k V. 



A* 



5 These are suppressed compared to the processes involving gluons by the multiplicity factor of 8 gluons 
versus 4 weak gauge bosons, as well as by a weak gauge coupling squared. 
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while the KK gauge bosons can decay into fermion pairs, — > Pi (always preserving KK 
parity). These processes are controlled by renormalizable interactions and the rates are very 
large compared to the Hubble scale at the times of interest. The exceptions to the above are 
i! which must decay into r' with a rate suppressed by A n and the radion KK states whose 
interactions are also controlled by A r . Of these r' is stable and by assumption remains out 
of thermal equilibrium. Since, as we will see, this assumption requires that A r be very large, 
the t' are quasi-stable, and freeze-out when their self-annihilation cross section becomes 
ineffective at xp = m t '/Tp. By this time, the number density of all the heavier KK particles 
is highly Boltzmann suppressed, and any r' from decays of the few remaining heavy particles 
give a negligible contribution. Therefore, we can integrate the Boltzmann equation including 
all processes until the time of t' freeze-out, assuming thermal equilibrium distributions. This 
will give a r' number density that we will call rffi. At a somewhat later time, the remaining 
if will decay out of equilibrium, giving an additional contribution to the r' number density. 
Assuming that there is no significant entropy production between the time of t' freeze-out 
and when the t' decay into r', we can effectively include this contribution at x = xf by an 
additive term rif) = n t '(xp). 

The total r' relic density today is then given by Q r > = m r in r i(xo)/ p c where p c = 1.05 x 
10~ 5 h 2 GeV cm" 3 = 8.06 x 10 -47 h 2 GeV 4 is the critical density of the universe today, 
x corresponds to the temperature of the universe today, and n r >(x ) can be connected to 
n r i{xp) = n$ + rffi from conservation of entropy. Thus, 

fi« = — ^T^TTTT 1 ^ [ XF dxx 4 C[x], (15) 



p c s(x F ) x F H(m t >) 



J-R 



where the entropy is given by 



45 x a 



and s = 2889.2 cm" 3 w 2.22 x 10" 38 GeV 3 is the entropy density today. Since T F ~ m t //30 
will be at most C(100 GeV), the relativistic degrees of freedom are those of the SM, e.g. 
g*{T = 100 GeV) = 92.25. 

The contribution coming from the decays of t' after freeze-out, as in the SWIMP sce- 
nario [HI H2], takes the form [see Eq. Q] 



(2) 2 m r , 1.04 x 10 9 GeV 1 x F 1 

m t > Mpi y/g^ {VCTt'i'^qq,gg) 

where vo t ti'^qq m is determined by the QCD interactions of t', as given in Eq. ^ (see also 
Fig. [5|. Neglecting all the quark masses (even the top quark mass gives a small correction 
for typical parameters), we have 

. 43 7ra., 

\VVt>i'->qq,gg) 



27 m 2 



where 43/27 = 7/27 + (2/9) AT/, with Nf = 6. The subleading corrections in the relative 
velocity, v, give also a negligible effect. The freeze-out temperature, parameterized by x F , 
is given by Eq. ^ with (vcr r v) - > ( vcr t'f->gg,gg) ; 9 — 2 x N c = 6, and m r > — > m t i. 
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Figure 8: Production of r' in the decays of (KK-parity odd) KK fermions and KK gauge 
bosons. Reversing the fermion arrows gives the production from anti-fermion decays. 



6.1.1 Decay Processes Before t' Freeze-out 



We now discuss the different contributions to the collision operator in Eq. (15), starting with 



the production of r' in the decays of KK-parity odd fermions and gauge bosons (see Fig. [8]). 
Using decays of KK quarks (which includes the top) for illustration, the contribution to the 
collision operator, per spin and color degree of freedom of the q\ is 



^Decay 



-E i/T -p 

e i J/ 1 



d 3 p, 



V 



q> — >r'q k 



(27T) 



where the decay rate is 
1 



2E. 



qJ 



(2Tc) 4 5(p q3 ~p r > -p k )\M, 



d 3 p r i d 3 p g k 
(27i) 3 2E r , (27i) 3 2E qk 



(19) 



(20) 



which is related to the decay rate in the CM via the time dilation factor 7 = E q j /rrij by 
rV-M-'gfc = Since the CM partial decay width is simply a constant, we can do the 

integral in Eq. ( 19 ) to obtain 



4?r m 3nCM ^ 

Air 



3 -pCM 



(27T) 



x 



-3/2 - 

2*> 6 



X 3 



X > 1 



< 1 



(21) 



(22) 



where K x {x) is a Bessel function, = nij/T = (mj/m t i) x, and now stands for any 

,, or r^jf G± according to the case (we give their 



„f -pCM -pCM -pCM -pCM 



expressions in Eqs. (71), (72), (76) and (77) of Appendix |BJ). We see that when Tr < rru 
(i.e. Xj > 1 for all relevant times), the contribution from C^ cay to the r.h.s. of Eq. (15) 
is exponentially small. If, on the other hand, Tr 3> mj, and taking into account that 



rrij > m t i ^> Tp, the relevant integral in Eq. (15) gives 



Xp 



dxx 4 ^^ 



e t',j 



CL3j j OC j 

mj/T R 



c5 



3tt 

Y 



(23) 



where e^j- = mt'/rrij, and therefore the contribution from decays is rather insensitive to the 
reheat temperature in this case. The support of the integral is 1/2 < Xj < 7, and we see 
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Figure 9: Tree-level scattering diagrams for gqi — > r'q k . Reversing the fermion arrows one 
obtains the diagrams for gq> — > r'q k . The contact diagram is required by gauge invariance. 

that the bulk of r' produced in decays of heavy particles occurs for temperatures of order 
the mass of the decaying particle, rrij. In particular, states with masses much heavier than 
the reheat temperature give only an exponentially small contribution. Note also that, since 
parametrically V™ ~ m 3 /A 2 , for states lighter than Tr we have f** dx x 4 Cj^ ay ~ m^m,/ 'A 2 ., 
and the individual contributions scale linearly with the KK mass, rrij. It follows that the 
contribution to the r' relic density from KK particle decays grows linearly with the reheat 
temperature, which sets the effective cutoff for how many KK states give a non-negligible 
effect. 



6.1.2 Scattering Processes 

We consider now the scattering processes. The contribution to the collision operator from 
each individual process gqi — > r'q k and for each spin/polarization/color degree of freedom in 
the initial state takes the form 

^scatt. _ f ~(E j+E g )/T l» 771 \ d 3 p g d 3 p q j 

C * = J 6 ( 4 ^-v) {2n)32Eg {2n)32Eqj , (24) 

where the scattering cross section is 

1 [ '/o \4r/ |2 d 3 p r > d 3 p q k 



AF j (27T) 6(p g +p q] - Pr> - Pqk )\M gg ^ r/qk \ {2n)32Eri , (25) 



and the flux factor is F = p qJ -p g (= V s IpI; with \p\ the three-momentum of one of the initial 
state particles in the CM). The relevant diagrams are shown in Fig. |9l and the analytic 



expression for the cross section is given in Eq. ( 78 ) of Appendix [Bj Following Ref . [43J , 
it is convenient to write d 3 p g d 3 p qJ = 27T 2 E g E qJ dE + dE_ds, where E± = E g ± E q j, and 
s = m 2 j + 2E g (E qJ — p qJ cos 9) is the total CM energy squared, while 6 is the angle between 
the momenta of the initial state particles. It is assumed that s > (m r i + mk) 2 - One finds that 
the physical region corresponds to E + > yfs, while \E- — E + m 2 /s\ < \/ E\ — s (1 — m 2 /s). 
The integrand is independent of E_, so that the E_ integral can be trivially done. The 
cross section (times the flux factor) depends only on s, and the integral over E + gives 
dE + ^ E\ — s e~ E+ / T = A /iTA' 1 ( v /s/T), where K\ is a modified Bessel function of the 
second kind. Thus, we can write 

r°° 87r 2 

x)Ctf*-[xi\ = m} y —KtMiu* - x 2 )(Fa gqj _ r , qk )du , (26) 



17 



where u = y/s/T and Xj = rrij/T are dimensionless variables. The full collision operator 
associated with scattering of KK fermions is then |^] 



G 



scatt. 



X\ 



AfC 



scatt. 



(27) 



where M = 2 (gluon pol.) x (A 2 — 1) (gluons) x 2 (q j spin) x N c (q j color) x 2 (g J and g J ) = 
8A c (iV 2 — 1), with N c = 3. For leptons and weak gauge bosons the multiplicities are different, 



but as we explained at the beginning of Section |6.1[ we neglect such contributions. We 
consider each SM quark flavor separately in the sum. Also, the allowed processes are of the 
form g + q~ — > r~q + or g + q + — > r~q~ , with ± referring to a KK-parity even/odd state. 

In the remainder of this subsection, we provide analytic approximations for the elemen- 
tary building blocks, C^ c fc att '. The reader mostly interested in the physics results may wish to 



jump to Section 6.2, where we discuss them. 



In the non-relativistic regime (T < rrij, which implies that the integral in Eq. (26) is 



dominated by the region u ~ Xj), we have the approximate result 



Fa 



non— rel 
gqi—tr'q* 



a. 



8N C A 2 ( u 



2 - x]y 



+ *? (r,/T) 2 



tan 



Ae 



1+Ae 



- 1 



x 9(Ae, e hJ ) [(G^) 2 + (G^) 2 ] (1 + Ae) + 4G™G™ e Kj 



(2? 



where 6*(Ae, ek, 



(l + Ae) 2 -4e 2 , Ae^e 2 



:,j ) — \l \^ 1 L -* K -I ""fcj' ^ c — c fcj c r' ' j) ^k,j — ""K/""]) ci±±u. y-r' ,j — ' ' "r' / ' • "J i 

i.e. we normalize with respect to the incident P mass. The Gf^, and Gf£ are the relevant 
couplings of the KK-radion (with % = 1 for r') to the fermion KK modes [see Eqs. (59), 



mk/rrij 



, and e r >j 



/rrij 



(60) and (54)]. This result is dominated by the s-channel diagram of Fig. |9| and we have 
explicitly included the qj width, Tj, which regulates the limit where the incoming gluon is 
ultrasoft and the s-channel q 1 is nearly on shell. 

In the relativistic regime (i.e. T 3> rrij) one can instead use 



Fa vel 



0£. 



rn- 



>r'q" 



64A C A 2 



[(G^y + (G^ 



r ijk) 



41og|£)-3 



(29) 



s ^> rrij of the exact cross section given in Eq. (78) of Appendix |B 



where we took the limit 

and defined Xk = m^jT . This result is dominated by the t-channel diagram in Fig. [9j The 
logarithmic divergence when — > is similar to the usual forward scattering singularity 
associated with the t-channel exchange of a massless particle. In Fig. 10, we show the 



exact scattering cross section as given in Eq. (78) of Appendix |Bj together with its non- 
relativistic and relativistic approximations, Eqs. (28) and (29), respectively. We choose 



typical parameters in order to illustrate the precision of the approximations. We will use the 

6 Although the couplings of the KK-radion to fermion zero-modes are suppressed by the corresponding 
fermion mass, its couplings to the KK states are similar for all the light family KK towers. Since the zero- 
mode contribution gives a negligible effect, we can describe the effect of the light families by a multiplicity 
factor, which would be appropriate in anarchic scenarios. 
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Figure 10: F x a(gP — y r't k ) as a function of z = m 2 /s, assuming m r i/rrij = 0.20, uik/mj = 
0.67, Gfff, = 0.50 and G±Q = 0.41, which correspond to the decay gt 5 — y r't 4 of the "small 
warping" benchmark scenario. 



The solid (blue) curve is the exact result given in Eq. (78) 
of Appendix the dashed (red) curve is the non-relativistic approximation of Eq. (28), 



and the dotted (green) curve is the relativistic approximation of Eq. (29). We also assumed 
Tj/rrij = 1/10, which is important only for z ~ 1. 



exact result in our numerical analysis, but the approximate limits provide a more transparent 
analytic understanding. 

For instance, in the non-relativistic limit, the u-dependence of the integrand in the colli- 



sion term of Eq. ( 26 ) is given by 

gmidu 



i^i(«)(« 2 



X 3 



I) 2 + Ij 



du 



(30) 



where jj = Tj/nij. It is instructive to rewrite ^nr as a function of 



w 



u 2 — x 2 



du 



2yjl + WJj 



dw , 



so that 



gNRdu 



1 



X 



;/m(.,' ;V m • „-,•) 



w 



dw 



(31) 



(32) 



This function determines the support of the scattering collision integral as follows. For 
Xj ^> 1 and writing u = Xj + Aw, we see that in this ultra non-relativistic regime (T <C m 3 -), 

jl. e -(%+A«)^ where the use of the asymptotic form 



we have K\{xjy/\ + wjj) = K\{u) r 

of the Bessel function is well justified. The exponential damping e~ Au ensures that the 
support of the integral satisfies Am < 0(1), and we can truncate the region of integration 

(and 

(33) 



at w m3X ~ 2/(^7j). Thus, using g^ndu 
otherwise), we can write 

« 27T 2 



5/2 _ 
2 X 3 e 



w dw 
w 2 +l 



for < W < Wr, 



x^C" 



non— rel r 



X.; 



(27T) 



7T 5/2 

2 Xj 



e Xj /nr (rrij , m k , m r > ) log 



1 + 
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where one should recall that Xj = (mj/m t ') x, and we defined 



fNR{mj,m k ,m r >) 



X : 



x Fa n ™- T fl k , 

gP—>-r't K ' 



(34) 



(u 2 - x 2 ) 2 + x1jj?_ 

a dimensionless function of the masses and couplings only (not temperature) that can be 



read from Eq. ( 28 ) . 



In the relativistic limit, x • <C 1, the -^-dependence of the integrand in the collision term 



of Eq. ( 26 ) is given by 



9r 



K x {u){u 2 -x 2 A x 



4 log 



u 



x k 



,„„. e « § + 1/log 



(35) 



. We may then 



This function has a maximum approximately at u 

replace u by u max inside the logarithm in g R . Using Ki{u){u 2 — x 2 ) = 2xjK 1 ( y Xj) ~ 2 for 
Xj -C 1, we can then write 



4 



16^ 



fR{mj,mk,m r >) 



4 log 



(36) 



where x-,- = (jrij/mt') x and = (mk/m t ') x, while 



64iV r . 



ljk) ' K^ljk) 



(37) 



is temperature-independent. Thus, the contribution to r' production via scattering of highly 
relativistic particles to the l.h.s. of Eq. (15) scales like x ~ 1/T (up to logarithms). Unlike 
the decay processes, scatterings that occur all the way up to the reheat temperature, Tr, 
contribute uniformly in the x-variable. These results follow essentially from the phase space 



for both cases (taking into account the time dilation factor in the decay widths). Eqs. (33) 



and (36) provide the basic analytic results to understand the production of r' via scattering 



of KK fermions off the plasma. 
In the left panel of Fig. 



11 



we show xlj C| c fc att ' [xj] as given in Eq. (26) for the process 
gP — > r't k , as a function of Xj = rrij/T, and factoring out the overall parametric dependence 

(27tW ft J? ( sc, hd, blue curve). We also show the non-relativistic (dashed, red) and relativistic 
(dotted, green) approximations given in Eqs. (33) and (36), respectively. We see that the 



relativistic approximation works well up to Xj ~ 0(1), while for larger Xj the non-relativistic 



approximation may be used. In the right panel of Fig. 11 we show the integral of the previous 



quantity with respect to Xj from Xj = tjij/Tr to Xj = rrij/Tp (solid, blue curve). Here Tp 
is the temperature when t' freezes-out, and it is assumed that mj/Tp ^> 1. Note that this 



integral differs from the integral appearing in Eq. (15) for the relic density by an overall 
factor of e\, • = (m t i /rrij) 5 . The dashed, red curve corresponds to an approximation where 



the simpler relativistic and non-relativistic expressions are used for Xj < x tr k and Xj > x tr 



respectively, for a conveniently chosen transition point xj k . We find that xj h = Srrij / (2m k ) 
gives a reasonably good approximation throughout the whole range, except in a region 
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Figure 11: Left panel: Collision operator for the process gP — > r't as a function of 
Xj = rrij/T. The solid (blue) curve is the exact result with the full cross section, Eq. (26), 
the dashed (red) curve is the non-relativistic approximation, Eq. (33), and the dotted (green) 
curve is the relativistic approximation, Eq. (36). Right panel: the integral w.r.t. xj of the 
previous quantity, as a function of mj/Tn (see discussion in the main text). It is assumed 
that mj/Tp ^> 1, where Tp is the t! freeze-out temperature. The solid (blue) curve cor- 
responds to the exact result. The dashed (red) curve is obtained by using the relativistic 
approximation when Xj < x 



x 



tr 
] k 



tr 

).k 



and the non-relativistic approximation when Xj > xj k , 



where 



2>mj / (2m k ) . The inset shows the same plot with a linear vertical scale. We have 



assumed the same parameters as in Fig. 10 



where rrij/T^ ~ few. In most of this region, however, the contribution to the final r' 
relic density is suppressed since P becomes non-relativistic. We also see from the rather 



weak logarithmic dependence on Xj of the relativistic expression given in Eq. (36) that, for 
rrij/Tji <C 1, the contribution to the integral in the l.h.s. of Eq. (15) becomes essentially 



independent of Tr [see also inset in the right panel of Fig. 11, and discussion after Eq. (37)]. 
This is similar to the decay process discussed in Subsection |6.1.1| We can therefore easily 



compare the contributions form scattering and decays for reheat temperatures that are large 



compared to the given particle masses. Using the result of Eq. (23) together with Eq. (21) 

see Eqs. dTTh or ([721) 



and 



m 3 J(167iA 



we have 



.(:;;<!■<■ -i- 1 (-try. 

E*dxx*Cfi«-[x] 



[47r/(27r) 3 ] [mf/(167rA2)] (3vr/2) 
?7r 2 /(27r) 6 ] (a 8 /N c ) (mf/A^) x 10" 1 



75tt 2 N c 
20a, 



(3? 



where the 10 1 is read from the right panel of Fig. 11 This estimate indicates that the 
contribution to Q r > from the decay processes dominates over the one due to scatterings. 



6.2 Constraints on a Superweakly Interacting KK-Radion 

In the previous subsection we have setup the formalism to compute the production of r' 
in decays of heavy particles, and in scatterings against the plasma. We established that 
the decay processes dominate over the scattering ones, essentially as a result of phase space 
considerations, and the fact that at any given energy there are KK states with masses 
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of that order, whose decays are important. Thus, we can focus on the production of r' 
via heavy KK mode decays. For a given reheat temperature, Tr, decays of particles suf- 
ficiently heavier than Tr give an exponentially small contribution, while decays of parti- 



cles lighter than Tr give a contribution to the integral in the l.h.s. of Eq. (15) of order 
(An) / (2n) 3 m 3 j TfJf(m t >/mj) 5 (3iT/2) = (3G 2 /647r 2 )mf/m : ,/A^ where G is a dimensionless ef- 
fective coupling constant for the corresponding decay vertex. The linear dependence on rrij 
translates into a linear dependence on the reheat temperature, which sets the effective cutoff 
for how many states give a non-negligible contribution to the r' relic density. Thus, the con- 
tribution from each KK tower is dominated by the states whose masses are around Tr. In 



fact, our discussion in Subsection |6 . 1 . 1 1 on the r' production due to decays of non-relativistic 
particles, shows that this contribution can be significant even for masses an order of magni- 
tude (or so) above Tr. One should also include a multiplicity factor for the internal degrees 
of freedom for the decaying particles. For instance, for each of the Nf = 6 quark towers, 
the multiplicity is 2 (q and q) x 2 (spin) x iV c (color). Thus, from the quark towers alone, 
setting rrij ~ Tr, and using Eq. (15), we can estimate the KK radion density arising from 
decays of KK modes as 



A, / V 10 TeV/ V300 GeV 



where we took g* = 86.25 Note that Eq. (39) is independent of m t i. A more precise de- 



termination is shown in Fig. [T2j were we show the region selected by WMAP in the Tr-Aj. 
plane, by treating the top towers separately from the light fermion towers (as suggested by 
the anarchic picture of flavor), and including the decays of KK gauge bosons, as well as the 
much smaller r' production by scatterings. The internal degrees of freedom associated with 
the top tower give a factor 2 x 2 x N c , while for the light generation towers we have a factor 
5x2x2xiV c +(3 + 3/2) x 2 x 2, which includes the five lighter quarks, three charged 
leptons and three LH neutrinos. The corresponding factor for KK gauge boson decays is 
3 (AT? — 1) + 3 x 4. We also include the r' produced in decays of t' after its freeze-out, which 
we previously denoted by Vt^h 2 , although this contribution is very small (see below). 

^From our discussion in Section [3j one might expect that if the reheat temperature is 
too large, the correct metric background would correspond to a blackhole solution, and the 
KK modes, including the radion tower cease to have meaning (in the 4D dual picture, the 
theory would be in a deconfined phase). Naively, we expect the critical temperature to be 



of order fc e g, and therefore we show in the left panel of Fig. 12 reheat temperatures up 
to k e g ~ 1 TeV. We also show reheat temperatures as low as Tr = 15 GeV (which is 
about the t' freeze-out temperature), where the KK-radion density is still small compared 
to its equilibrium distribution. The radion decay constant that reproduces the observed DM 
density ranges from 10 4 — 10 16 GeV, depending on Tr. 

However, given that we do not have the corresponding blackhole solution that would 



allow us to compute the critical temperature, we also show in the right panel of Fig. 12 the 
result of assuming that the KK theory is the correct description up to temperatures of order 
its 4D cutoff (which we estimate to be of order 5 — 10 TeV). Even in this case, WMAP selects 
a A r not much larger than 10 16 GeV. For reheat temperatures above a TeV, lower values 
of A r would lead to overproduction of r', while larger values would make the r' species a 
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Figure 12: The WMAP constraint in the Ar-Tr plane. The left panel shows the allowed region 
on a logarithmic scale for 15 GeV < Tr < k c s ~ 1 TeV and 10 4 GeV < A r < 10 16 GeV. 
In the right panel we show, on a linear scale, the WMAP allowed region up to Tr = 10 TeV 
(the cutoff of the KK theory), and illustrate the effect of the number of KK modes included. 
The blue (upper) region includes contributions from KK states up to an order of magnitude 
above the maximum reheat temperature shown. The red (lower) region shows the result of 
truncating the KK tower at 10 TeV. 



subdominant component of the observed DM density. However, note that, within the model, 
A r cannot be arbitrarily large, but is bounded from above by the Planck mass. Thus, 
quite aside from the likely difficulty in nucleating to the true vacuum from a "deconfined" 
one if the temperature is above the critical temperature, this result suggests that Tr cannot 
be much higher than the cutoff of the effective theory, or else the KK-radions would be 
overproduced. 

As mentioned above, states somewhat heavier than Tr in general cannot be neglected. 
This raises the question of how to treat these states when Tr is close to the cutoff of the 
effective theory, even assuming that the theory remains in the "confined phase", so that a 
KK radion can be properly identified. Nevertheless, even if the relevant degrees of freedom 
above the cutoff are not accurately described by the KK theory, one might expect that 
their couplings to the KK-radion are still suppressed by A r . One can therefore get an idea 
by including the contributions from KK states up to about one order of magnitude above 



the cutoff. The result corresponds to the upper, blue region in the right panel of Fig. [12 
which exhibits the linear dependence on Tr previously discussed. However, we also show 
the result of truncating the KK tower at a cutoff of 10 TeV (lower, red region), as might 
be appropriate if for some unknown reason the relevant degrees of freedom above that scale 
are significantly more weakly coupled to the KK-radion than expected from the KK picture. 
The comparison between the two curves gives an idea of the effect of the heavy states at 
larger reheat temperatures. 



7 Recall that these scales are only effective scales, characterizing the strength of radion and graviton 
interactions. In the present "small warping" scenario, the physical cutoff of the theory on the UV brane is 
in fact significantly smaller, of order the 5D Planck mass, M5. 
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Figure 13: Left panel: t! lifetime as a function of A r for m t > = 600 GeV, m r i = 300 GeV 
and G = 1 (see text). Right panel: t' yield, Y t /(xF), as a function of m t i. For reference, in 
the right axis we show Vt t ih 2 = (m t i / m r /)Vffi h 2 corresponding to the yield shown in the left 
axis and the mass m t i given by the curve. 



In this scenario, the large value of A r / would make the NLKP long-lived. For instance, 
if the NLKP is t', and assuming m t i > m r i + m t , it decays via t! — > r't with a decay width 
given by 



G 2 m % 
16 - 



ITLf m? 



m t 



m t > 



(40) 



where G is an effective coupling constant [see Eq. (71)]. In Eqs. (76) and (77), we give 



the decay widths for the KK-parity odd Higgses in the case that one of these is the NLKP. 
The constraints are similar to the case of a t' NLKP, and in the following we concentrate 
on the latter case. In the left panel of Fig. 13, we show the t! lifetime as a function of 
A r for m t > = 600 GeV, and taking G — 1. The result for other values can be roughly 
obtained by simple rescaling, according to the prefactor in Eq. (40). These late time decays 



are potentially dangerous. When massive particles decay into high energy quarks or gluons, 
the latter rapidly fragment into hadrons, and if the decays of the mother particle occurs at 
times 10~ 2 s < t < 10 4 sj^jthe scattering off background protons and neutrons via the strong 
interactions can lead to a dangerous interconversion of background protons and neutrons 
even after the freeze-out time of the neutron-proton ratio. The effect tends to increase n/p, 
thus increasing the 4 He abundance with respect to the standard predictions of BBN. For 
lifetimes smaller than about 10 4 s, EM energy deposition is weakly constrained |44J . hence 



the left panel of Fig. [13] suggests that the constraints arise primarily from the hadronic energy 
injection. 



In the right panel of Fig. 13 we show the t! yield, Y t i(xp) = n t i(xp) / 's(xp), which is set 



by the QCD processes of Eq. (18), and depends only on m t i. In the right axis of the same 
figure we show the associated Q t 'h 2 , which leads to Q^h 2 = (m r / /m t >)Qt>h 2 , according to 



8 At earlier times, such hadrons predominantly scatter off the background photons and electrons, trans- 
ferring most of their kinetic energy and reaching kinetic equilibrium. 
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Eq. (17). We can see that the latter one is typically of order 10~ 2 or less, and therefore 
much smaller than the total Q r rh 2 . This suppression in the t' yield allows the picture to be 
roughly consistent with BBN constraints even when the t' decay during or somewhat after 
BBN. The yield can be constrained as a function of the lifetime of the long-lived relic and its 
mass, which determine the amount of injected energy and when this injection occurs [4"5ll4"6] . 
We can therefore bound A r as a function of m t > as follows. For a given my the right panel 

151 we determine how late the relic 



of Fig. 13 gives the yield, Y t i. From the results of Ref. 
can decay without being in conflict with BBN constraints (typically around 200 s). This 
can then place an upper bound on A r via the t' decay width. Assuming for illustration that 
G 



1 and my — 300 GeV in Eq. (40), we obtain in this way 



m t > [GeV] 
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which gives an upper bound on A r for a few typical values of my (or, equivalently, a lower 
bound on my for several values of A r ). Of course, if my < m r / + m t , so that we have a 
three-body decay, the BBN constraint becomes more stringent. 

There is a stronger constraint from the late decays of the radions (produced in decays of 
KK states, analogously to the production of the KK-radions described above). In fact, since 
the masses and wavefunctions of the radion and KK-radion are nearly identical, their cou- 
plings to KK modes are very similar (with appropriate replacements of the semi-degenerate 
even and odd modes, as dictated by KK parity). We have checked that the number density 
of radions is very similar to that of KK-radions. The radion yield can then be expressed as 



SQlTL r i 



(41) 



where p c and so are the critical density and entropy of the universe today, respectively, and 
Q r i is the DM relic density. For instance, for m r > = 500 GeV, and assuming Q r /h 2 = 0.1, we 
get a radion yield of V r ~ 7 x 10 -13 . On the other hand, the decays of the radion into W , 
Z, and h + give a radion lifetime of 



ml 



-i 



0.1 x 



A,. 



10 15 GeV 



500 GeV\ 3 

I sec 

r i J 



m r 



(42) 



Here we assumed that mw,z, m h m r and neglected the smaller width into tt and massless 



reference parameters (m r 



gauge bosons (we give the exact formulas in Eqs. (84)-(87) of Appendix B.7). For the above 

~ 13 ), Ref. [15] gives an upper 

we 



500 GeV with Y r « 7 x 10' 
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bound on the radion lifetime of about 0.5 sec. Referring to the left panel in Fig. 
infer that the reheat temperature should be at most in the few hundred GeV range, but 
that under this assumption it is possible for KK-radions to account for the DM relic density, 
while being consistent with BBN constraints (from both radion and t' decays 0). Recall 

9 Heavier KK radions are more weakly coupled and therefore decay even later. However, they are also 
produced in smaller quantities, and given that the reheat temperature is significantly below their mass, we 
expect their number density to be exponentially suppressed. 
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that a maximum reheat temperature of order a TeV may be required in order to not exceed 
the critical temperature for a deconfmement /confinement phase transition, as discussed in 
Section [3j 



7 Direct and Indirect Detection 

We now turn to the feasibility of KK-radion detection in either direct or indirect DM searches, 
as well as collider experiments. 

Direct detection: 

As has been mentioned before, most of the interactions between the KK-radion and the SM 
particles are non-renormalizable, and suppressed by A r . After EWSB, renormalizable interac- 
tions with the SM can be induced by terms of the form \ [5(y — L) + 5(y + L)) ^gind CH^H^H 
(subject to the KK-parity symmetry). For simplicity, we restrict here to an exactly IR lo- 
calized Higgs field, so that the previous operators are IR localized and involve the Ricci 
scalar constructed from the induced metric. These induce kinetic mixing between the radion 
and the SM-like Higgs [MJ |33]. In the present context they also induce kinetic mixing be- 
tween the KK-radion and the KK-parity odd (CP-even) h_ that resides in the inert Higgs 
doublet, as already mentioned at the end of Section |4j Canonical normalization is achieved 
by the field redefinitions r± = r±/Z and h± = h± + (5/Z)f±, where S = 6y/2£,v/A r and 
Z 2 = 1 + 12£(1 — 6£,)v/A r , with v = 174 GeV. In the following, we will use the notation 
r + = r and r_ = r' to denote the radion and KK radion states. After the previous field 
redefinitions, the mass matrix becomes non-diagonal. The mass eigenbasis is obtained by 
the orthogonal transformations: 

h± ) \ U h±jL U h±iH J V <t>±,H 

where the subscript L (H) refers to the lighter (heavier) mass eigenstate within the KK- 
parity even and odd sectors, which we denote with the + and — subscripts, as usual. The 
matrix elements above are explicitly given by 

a± - c± - A± a±-c± + A± 

U r±,L — = = ll0 i U r± ,H 



(43) 



VW - c ± - A ± )2 + 46| ' ±! v /(a ± -c ± + A ± )2 + 4fe2 



2b± 2b± 

Uh ± ,L = ,. = tl0 , Uh ± ,H = I, = = Al0 i (44) 

V(a± - c± - A±) 2 + Ab 2 ± v /(a±-c± + A ± ) 2 + 46| 



where A± = y/(a± - c±) 2 + 46^, and N 



a± = ^ ( m r± + 62m l ± ) > b± = £ ™,h± > c ± = m l ± ■ (45) 

Here m 2 ± and m 2 h± are the mass parameters in the (r±,h±) basis. The physical masses in 
the KK-parity even and odd sectors are given by 

m ±,L = \ (a± + c ± - A ±) , m ±,H = o ( a± + c± + A± ) • ( 46 ) 



10 Note that, when a± < c ± , {U r±tL ,U h± , L ) -> (—1,0) and (U r± , H , U h± , H ) -> (0,1) as 6± -> 0. When 
a± > c±, one has instead (U r±j z, Uh ±i £) — > (0,1) and {U r±y H,Uh ± ,H) — > (1,0) as 6± — > 0. 
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The DM particle is the lightest KK-parity odd state above, 0dm = 4>-,l, with mass 
hidm = rri-,L- Now we can identify the interactions relevant for scattering of the DM 
candidate against nuclei. As shown in [IB] , there are r_h + h_ and r + h_h_ terms of the 
form 



A -r_ [dnh+Wh- - Am 2 h+ h + hJ\ ->• ^ — U r _ >L U h+t0l U h _ >L (D0 +iQ + 8m2 + +>a ) 



DM > 



2 



[d^h-d^h- - 2ml_h-hJ\ -»■ 

" E !T U r + ^UL,L [n0+,a + (2m DM - 4mL)0+,a] 



a=L,_ff r 



" I VDM 5 



where, on the r.h.s, we isolated the terms quadratic in 0dm- To obtain the above expressions 
we integrated by parts and assumed that 0dm is on-shell: D0dm = — m DM^DM- There is 
also a contribution from the Higgs potential 



L,L l C+,«TOM 



The exchanged KK-parity even state can be either +i l or 0+,//. The corresponding Feynman 
rules for the +i l0 dm and +i #0 DM vertices are, respectively, (2i/ A r )(ZL } n q 2 +M. 2 L H ), where 
g is the 4-momentum of the KK-parity even state, Z a = U r _,LUh + , a Uh_,L — U r+ta Ul_ L , and 
Ml = &ml + U r _ iL U h+ , a U h _ iL - (2m DM - 4mLR+,a^L,L - (3/2V2)(A r /v)m 2 h+ U h+ ' a U 2 _ iL . 
Identifying also the Higgs (i.e. h + ) component in + i and 0+ ! h, we obtain the DM-nucleon 
scattering cross-section 

<W^ DM * - |f , , ™* . (U h+>L Ml m l H + U h+ , H MW + , L fm 

nA r m +L m +H {m N + mY)M) 

In the above, we neglected the coupling of the radion to the nucleon compared to the Higgs- 
nucleon coupling, g hNN rs 340 MeV/246 GeV [37]. Note, however, that for A r ~ 2 TeV, the 
radion-nucleon coupling, of order m N /A r , need not be completely negligible (it can be easily 
incorporated if wanted). 



In the left panel of Fig. 14 , we show contours of constant c</> DM Ar->^ DM JV in the £-m r _ plane, 
taking A r = 2 TeV and rrih_ = 1 TeV. Note that, for practical purposes, the DM physical 
mass is very close to m r > = m r _ for m r i < m h _, and very close to rrih_ otherwise. Thus, the 
horizontal axis in the plot can be readily interpreted in terms of moM- The dashed (white) 
lines correspond to constant \U r _^\ 2 (the r_ content of the DM matter). As discussed at the 
end of Section |4j the h- content, given by 1 — \U r _ t i\ 2 should be small, or else the DM relic 
density will be too small due to efficient annihilation via the Higgs couplings. However, we 
see that there is a region where the DM is mostly r_, with a sizable nucleon cross section 
(e.g. below the dashed line of \U r _^\ 2 = 0.9). In fact, the current direct DM searches may 
have a sensitivity to this region, as shown by the gray region corresponding to the most 
recent exclusion by the XENON100 experiment 
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Figure 14: Left panel: Contours of constant 0> DM jv->0 DM iv m the £^m r _ plane, as given in 

2 TeV. Also shown is the region currently excluded by XEONO100 

2 (with 



Eq. gTp, for A r 

(gray). The dashed, white lines represent the r_ content of the DM candidate, \U r _ t L 

12 _ 1 |7T 12 



Uh-,L\ 2 — 1 — \U r _ ) L\ 2 being the h_ content). Right panel: Nucleon cross section as a 
function of m r _, for different values of A r , and for constant £ 
to rrih_ = 1 TeV. 



0.1. Both figures correspond 



In the right panel of Fig. 14 , we show curves of the nuclear cross section for several values 



of A r for fixed £ = 0.1 and rrih_ = 1 TeV. The plateau for m T > > rrih_ corresponds to the 



limit in which the DM is pure h_, where Eq. (47) reduces to 



^umN^umN ~ -73— T {mN + mDM )2 . I 48 ) 

with mDM = rrih_ = 1 TeV (we also assumed A r ^> v, so that the exchanged h + has no 
radion admixture). 

Indirect detection: 

Regarding indirect detection, the most interesting annihilation channels are those that lead 
to positrons, neutrinos or photons. The direct annihilation cross section into e + e~ pairs 
takes the form vcv r /_> e + e - ~ (Gf 100 /47r) (rn^/Af), where Gnoo is the coupling for the contact 
interaction between two KK-radions and e + e~ (see second diagram in Fig. [I]). In anarchic 
scenarios, due to EWSB effects and the subsequent mixing with the electron KK resonances, 
the coupling Gnoo can be as large as O(30). However, the ml factor makes this channel 
extremely suppressed. Mixing with the Higgs does not help since it also couples to elec- 
trons via the small electron Yukawa coupling. Similarly, the direct annihilation rate into 
neutrinos is negligible. Softer positrons/neutrinos can be produced in the decay products 
of the dominant tt and Higgs (or gauge boson) annihilation channels, but these the fluxes 
are expected to be small. For instance, given that the spin-independent DM-nucleon cross 
section is bounded from above by about 10~ 43 cm 2 ps 10~ 7 pb for mDM ~ 300 GeV, neutri- 
nos from DM annihilations in the Sun would produce at most a couple of neutrino events 
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Figure 15: Tree-level (from heavy gauge boson exchange) and 1-loop (from KK fermions) 
contributions to r'r' — > 77. The ± indicate the allowed KK-parity combinations. 



per year [IH] at neutrino telescopes such as IceCube/Deep Core. This is well below the 
atmospheric neutrino background, and the detection in this channel is not likely. 

We therefore focus on the photon flux from KK-radion annihilations, starting with a 
possible line signal. There are tree-level direct annihilations from the exchange of KK- 



photons, as shown in Fig. 15 A given KK-parity odd photon induces the effective dimension- 
8 operator —2g\- / [Is^m 2 ) d^[r' F^ v ]d a [r' F av ], where gij is the coupling between the KK- 
radion, the KK-photon with mass rrij, and the massless photon, as given in Ref. [IB]. Since 
both A r and nij may be expected to be of order a couple TeV, this contribution is small 



(we give the full cross-section for any mj in Eq. (68) of Appendix |BJ). There are also 1-loop 
contributions, e.g. from fermion KK loops, as shown in the right diagram of Fig. 15 Imposing 
only [7(1 ) em gauge invariance, the effect from each KK level is logarithmically divergent by 
power counting, which would seem to lead to a linear divergence after summing over KK 
modes, rj Nevertheless, as pointed out in Ref. [50] for the case of the linear radion couplings 
to 77 (where a triangle rather than a square diagram is involved), these divergences simply 
renormalize the 5D operator ^Jg~ FmnF mn . Unfortunately, this operator does not lead to 
interactions involving two KK-radions and two gauge zero-modes, as pointed out in Ref. [IB] • 
The remaining finite 1-loop pieces contribute to the dimension-8 operator quoted above, and 
result in a negligible effect. 

Nevertheless, one can still write 5D operators that do lead to the (r') 2 77 interaction if the 
stabilizing scalar field is involved (since the KK-radion is really an admixture of metric and 
stabilizing scalar fluctuations). However, we find that a more important contribution can 
arise from an IR localized operator involving instead the extrinsic curvature. For instance, 
consider the operator [5(y + L) + 5(y — L ) ] ^ ind (i]/8Al)K 2 F^ V F^ U , where A 5 is the 5D 
cutoff, the indices are contracted with the induced metric, and K = 4[A'(y) + d y F(x, y)} is 
the trace of the extrinsic curvature, including radion fluctuations about the background, as 
parameterized in Refs. [501 EE] • Reduction to 4D with the flat photon wavefunction, 1 / \/2L, 
as well as the KK-radion wavefunction, which satisfies F (L) m 1 and Fq(L) w 2A'(L), leads 



11 Although there is a double sum in the square diagram, the "diagonal" couplings, i.e. those for the pairs 
of KK-parity even and odd fermions with closer masses are largest, while the corresponding couplings rapidly 
decrease as the KK-fermions become more separated in mass 1 1 8 1 . 
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to the dimension-6 operator 



12 



e 2 K /2 



8A 2 



r l2 F^F^ , (49) 



where e 2 K ~ 64?7/(/c eff L) x (A; eff /A 5 ) 3 , with fc eff = A'(L). NDA suggests that the unknown 
dimensionless coefficient, r], can be as large as 37r/2 [SB]. Taking also A; e gL ~ 60 and 
A5 ~ 5A; e fj, one estimates e 2 n w 0.04. Notice that the same operator contains a term 
without radion insertions, which renormalizes the electric charge. However, this effect is 
very small and we neglect it in the following. 



The interaction of Eq. (49) leads to a direct 77 annihilation cross-section given by 



(a 27 v/c) « — —2 , (50) 

\ ^V r / m DM 

where ttidm = wv, and a = e 2 /47r is the fine structure constant. In the following, we do not 
include the additional 7Z channel, nor any mixing with the Higgs. 

The differential photon flux from the galactic center can be written as 

d ^< - 5.66 x 10" 12 cm-V 1 ^ ( ( ±™) ' J(AO)Afl , (51) 



dE dE \ 1 pb J \ niDM 

where dN^/dE = 2 x 5(E — m-^yi) is the differential photon yield in the 27 annihilation 
channel, J (An) = (1/Afi) J An J(ijj)dfl integrates over the angular acceptance of the detector 
AQ, and J (iff) is conventionally defined as 



JW '> = a5kp^( o.3GeV/c m3 J J, «™ • (52 > 

where r 2 = l 2 + r 2 — 2/ro cosip, with ro ~ 8.5 kpc the distance from the Earth to the galactic 
center. The integration is along the line of sight, dl, and encodes the information about 
the DM distribution, assuming a spherical DM halo of energy density p(r). The quantity 
J (Aft) can vary over several orders of magnitude depending on the halo model. We will use 
here the Navarro- Frenk- White (NFW) profile [5T] . but the results for other profiles can be 
obtained by simple rescalingj^ We have also used the standard DM density at the Sun's 
location of 0.3 GeV/cm 3 , although it has been recently argued that it should be somewhat 
larger [S3]. We take AQ = 10~ 5 sr, as would be relevant for the HESS and Fermi-LAT 7-ray 
experiments. We then have J(AQ)AQ f» 0.15. 

We take into account the detector energy resolution by convoluting 



dE 



measured 



dE G(E, E )—j-j^- , (53) 



12 We point out that, in spite of being localized on the IR boundaries, the above operator leads 
to a non- vanishing linear coupling between the radion field and two photons: fd 5 x y/g F(x,y)T — > 
— Jd^x (877 k^ s /LAl)(r/A r )F^ u F^, where T is the trace of the stress-energy tensor induced by the op- 
erator. This result holds even if the gauge fields are strictly IR localized, and therefore can be interpreted 
as an explicit breaking of conformal invariance in a dual 4D picture. 

13 It has been claimed that cored profiles are the rule [55], in which case the photon flux would be signifi- 
cantly smaller. 
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Figure 16: Continuum photon flux and 2^-line signal from KK-radion annihilation at the 
center of the galaxy. For the line signal, we show three detector energy resolutions: 10% 
(short dashed, black line), 5% (solid, blue line) and 0.5% (dashed, red line). 



where, following [M], we take a Gaussian kernel, G{E,E ) = l/(V^E a) e -^-Bo) 2 /(2^) 
with a = £/2.3, and £ is the detector's relative energy resolution. We show the line signal in 



Fig. 16 for mDM = 300 GeV and A r = 3.5 TeV (the rest of the parameters are those for the 
"strong warping" benchmark scenario) . We also show in the figure the continuum 7-ray flux 
spectrum as computed by micrOMEGAs 2.4 [55], which also reproduces our result for the DM 
relic density (and, for the previous parameters, gives Q r rh 2 ~ 0.11). We see that the line 
feature is distinguishable from the continuum spectrum signal, perhaps even for the current 
typical detector energy resolutions of about 10%. 

It is interesting to compare to UED scenarios [56]. The continuum photon signal was 
computed in Ref. [57], where it was pointed out that the resulting spectrum is rather flat 
up to the DM mass. This is due to the large branching fraction for direct annihilation of 
the KK into lepton pairs, which results in a hard photon spectrum. In our case, with 
annihilations mostly into top pairs, gauge bosons and the Higgs, the resulting spectrum is 
softer than in the UED scenario. The photon line signal, which arises from a calculable 
1-loop effect, was also studied in UEDs [581 El]- It was found that the line peak could be 
easily seen provided the detector energy resolution was small enough. In the KK-radion 
scenario, we point out that an incalculable 5D interaction can also lead to a prominent line 
signal. One should also remark that the KK-radion can be lighter than which can be 
very interesting from the point of view of FERMI-LAT, which has sensitivity for photon 
energies between about 20 and 300 GeV. We find that the integrated photon flux from the 
KK-radion signal, in the above range, is $ 7 ~ 3 x 10~ 12 cm -2 s -2 , which is somewhat more 
than an order of magnitude smaller than FERMI-LAT's sensitivity of about 10 -10 cm -2 s~ 2 . 
However, one should recall that there are large uncertainties associated with the DM halo 
model profile that can change the signal by orders of magnitude. 

Collider signals: 

Finally, we comment briefly on the collider signals of our setup. Regarding the more generic 
KK phenomenology, which can be expected to share features of the widely studied RS 
scenarios (see [59] for a review), we only note that the generic 1st KK- level states (which 
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are KK-parity odd) and 2nd KK level states (which are KK-parity even) can lead to signals 
both with and without missing energy. In this respect, the situation is rather different from 
supersymmetric scenarios with R-parity. In addition, we have emphasized that KK-parity 
even and odd states come in relatively degenerate pairs, which allows to access both types 
of states almost simultaneously. This can be contrasted to the UED case with KK-Parity: 
although the first KK-level UED particles always lead to missing energy signals [60J , and the 
second level KK states can lead to pure SM particle final states [601 EE] , the latter are about 
a factor of 2 heavier than the former in 5D constructions [2] (in 6D constructions such as the 
"chiral square" of Ref. [I] this factor is just \/2, which can be further diluted due to radiative 
effects [62]). Note also that there are significantly lighter KK-parity odd fields (discussed 
next) that mitigate the requirement of pair-production, so that the reach in KK-parity even 
and odd states should be similar. 

Our main focus here is on the DM-related collider signals. Of particular importance are 
the KK-parity odd states such as the first excited top resonance (probably the SU (2) singlet), 
or perhaps the "inert" Higgs doublet. These states may be expected to be parametrically 
lighter than the remaining 1st KK-level states, and thus provide the most promising portal 
for r' production. For instance, in the "strong benchmark scenario" discussed in previous 
sections, we have m r i ~ 300 GeV and m t i ~ 600 GeV. The vector-like quarks have a strong 
production cross section at the LHC of a couple picobarn [63, 64j. Their main decay is 
i! — > tr', thus leading to a ti +$t signal. This topology was studied in a model-independent 
way in [65]. Given the large $t, which can be used for tagging, one can consider the case 
where both top quarks decay hadronically, thus allowing for their full reconstruction. It was 
shown in the previous study that a judicious set of cuts and kinematic variables can allow 
for a t' discovery at the LHC with less than 10 fb _1 , and in fact that both rrif and m r > 
can be statistically measured. Interestingly, within the context of our model, we can have 
access to the LKP mass by studying the closely related radion mode, that behaves in many 
ways like a Higgs. For m r ~ m r i ~ 300 GeV, the radion will be produced by gluon fusion 
and has a significant branching fraction into ZZ p4~l 135] 166] , thus allowing for a clean mass 
measurement in the fully charged 41 channel. We can therefore get an interesting hint for 
the expected degeneracy between r and r', and the relatively precise knowledge of the radion 
mass can be used to study in more detail the properties of the missing energy signal. 

It would also be extremely interesting to measure the radion decay constant, which sets 
the total width for both t' and the radion. For the above masses, and for A r ~ 3 — 4 TeV, 
as suggested by the relic density scenario discussed in Section |4| these widths are of order 
a fraction of a GeV. This may be too small to be measured at a hadron collider, but would 
be accessible in a lepton machine. Nevertheless, one might be able to infer A r from a 
measurement of the radion production cross section. This could be a first test of the required 
properties to account for the DM relic density. Eventually, we expect that a compelling case 
for the identity of DM could be made based on collider measurements. 

In previous sections, we also considered scenarios with radion decay constants much larger 
than the TeV scale. In the coannihilation scenario of Section [5j with A r as large a 10 4 TeV, 
the t' still decays promptly but the produced tops (and even the W's) are expected to be 
off-shell. The analysis of this signal is expected to be challenging. It may also be possible 
for the t' to decay, via intergenerational mixing, into the lighter families. On the other hand, 
the non-thermal scenario studied in Section [6] can potentially lead to a t' that is long-lived 
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on collider time scales, if A r > 10 10 GeV. The same applies to the "inert" Higgs doublet 
(which has a charged component) if it turns out to be the NLKP. These stable, charged 
states would lead to spectacular signals at the LHC. 

8 Conclusions 

In this work we have considered warped extra-dimensional scenarios with a KK-parity sym- 
metry that makes the lightest KK-parity odd particle stable and hence a DM candidate. 
The LKP is expected to be the first KK excitation of the radion field, which is assumed to 
be stabilized at tree-level by the Goldberger-Wise mechanism. The radion mass is expected 
to be of order a few hundred GeV, but its interactions -which are controlled by the radion 
decay constant A r - can vary over orders of magnitude. We considered a number of scenarios. 

Thermal KK dark matter 

1. The desired KK-radion relic abundance can naturally arise when A r is in the multi-TeV 
range. EWSB plays an important role here in that it induces strong mixing in the top 
KK tower, and as a result a relatively large effective coupling between the KK-radions 
and a tt pair that partially offsets the suppression from A r . 

2. The presence of light fermion resonances (e.g. the first KK excitation of the RH top) 
can, via coannihilation effects, deplete the number density of thermal KK radions down 
to acceptable levels for even larger radion decay constants of order 10 4 GeV. Provided 
the required degree of degeneracy is present (of order 10-15%), the relevant annihilation 
cross section may be purely controlled by the QCD interactions. 

3. Alternatively, a small KK-parity odd Higgs component of the DM candidate, can be 
rather efficient in enhancing the DM self-annihilation cross section. This is due to the 
relatively strong interactions of the CP-even (but KK-parity odd) KK Higgs with the 
fermion KK states. If this Higgs component is sizable, one expects that these KK 
particles would constitute only a fraction of the observed DM density. However, it is 
possible that a combination of radion-Higgs mixing and a large radion decay constant 
can result in the required annihilation cross-section to obtain the DM relic abundance 
via the freeze-out paradigm. 

Non-thermal KK dark matter 

4. In UED-like scenarios, where the bulk curvature is small compared to the 5D Planck 
mass, the radion decay constant can be as large as the Planck scale. If the KK- 
radion interactions are so weak that the KK-radions never reach thermal equilibrium, 
a lower bound on the radion decay constant can be set from the requirement that 
these non-thermal KK radions do not overdose the universe, for a given reheat tem- 
perature. In such a case, there is a population of long-lived radions and NLKPs 
that can decay around the time of Big Bang nucleosynthesis. Consistency with BBN 
puts an upper bound on A r , which controls their lifetimes. This constraints implies 
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A r < O(10 15 GeV), while the reheat temperature should be around the EW scale 
or somewhat below. We also noted that there may be another upper bound on the 
reheat temperature, of order TeV, to avoid a dangerous high-temperature deconfine- 
ment/confinement phase transition. The point is that in the scenario with very large 
radion decay constants, the rate of bubble nucleation of low-temperature vacua would 
be extremely suppressed, and furthermore the associated reheat temperature after the 
transition may not be high enough to allow for a successful BBN. These constraints 
can be avoided if the universe never reheated above the critical temperature, which can 
nevertheless be high enough to allow the non-thermal DM KK-radion picture above to 
be realized. 

We also explored the prospects for direct and indirect detection of KK-radion DM. If the 
scenario has a relatively low radion decay constant, direct detection is feasible and a region 
of parameter space is already excluded by XENON100. A small KK-Higgs component is 
essential. For indirect detection experiments, annihilation into positrons or neutrinos is 
expected to give rather small fluxes (both from direct annihilations, as well as in decay 
chains of the dominant ti /gauge /Higgs channels). However, there may be an interesting 
photon signal that includes a prominent photon line. This line signal can naturally fall in 
the energy range of sensitivity of Fermi/LAT, and the total flux may be measurable for cuspy 
DM halo models. 

We also point out that warped scenarios with KK-parity would lead to interesting collider 
signals, both with and without missing energy. In particular, one may be able to probe the 
characteristic degeneracy between the DM candidate and the radion mode (which leads to 
signals similar to the Higgs), at least if the first KK-parity odd excitation of the top is in 
the expected few hundred GeV range. The scenarios with very large radion decay constant 
can lead to highly ionizing charged tracks. 
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A Radion- fermion interactions and EWSB 

The most general interaction Lagrangian of a single and two fermions can be written as 

£*H, = y \9f k ¥ L inl + 9^¥ R iH\-(™f k VM + h.c.) , (54) 

ijk=0 r ^ ' 

_ *r> _ _ 

where tp fix = \ [^l^d^x ~ {.9^ip)l^x\ ? the and g^j} are dimensionless, and the have 
dimensions of mass (these coefficients satisfy the symmetry properties = g^, g^ = gf^ 
and m\R = mf^). The sum runs over all fermionic mass eigenstates (as well as the rj). In 
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general, the ipi are not gauge eigenstates, since the KK towers associated with different bulk 
fields mix when the Higgs gets a VEV, thus inducing EWSB mixing mass terms of the form 



EWSB 



E (™% L xl R X k 2 , L + h.c.) 



(55) 



jk=0 



Here \i an d X2 stand for two different bulk fields in the gauge eigenbasis (the SU(2)l 
singlets and doublets of the 5D SM), while the mj| L are proportional to the Higgs VEV, 
v. For a Higgs localized on the IR brane, these mixing mass terms are given by = 
v {Y hD /L)\ [1 + {-iy +k ] f{ R {L)fl L {L), where v = 174 GeV, and f^ L are the properly 
normalized fermion wavefunctions evaluated on the IR boundary at y — L, and the factor 
in square brackets reflects the KK parity symmetry of our setup. Perturbatively, the 5D 
Yukawa coupling, Y 5 d, is related to the 4D Yukawa coupling by Y±d = (Y 5 d / L) fl(L) f R (L) , 
where m^o = Y^v is the corresponding SM fermion mass. However, one should match to the 
observed fermion mass after diagonalization of the full KK mass matrix. This is especially 
true in "anarchic" scenarios where the 5D Yukawa couplings are large and mixing with the 
KK states can be important. 

There are linear interactions of the vi that have the same origin as the EWSB mass terms 



of Eq. ( 55 ) and that take the form 



/-EWSB 



ijk=0 



Ar 



^JfeXipX^ + h.c. 



(56) 



For Higgs doublets localized on the IR branes, one has X RL = 2(Y 5D /L) [l + (— x 
Fi(L) f( R (L) f2 L (L) , where Fj is the rj wavefunction. The coefficients vX^, after rotation 
to the mass eigenbasis, can be assumed to be part of the of Eq. (54). The remaining 
contributions to the coefficients in Eq. (54) arise from the bulk terms inuie 5D action, and 
are obtained from the overlap integrals given in Ref. |18j . 14 after properly applying the 



unitary transformations that diagonalize the mass matrix, + fh RL = Ur mp^ysUj , where 
m KK contains the KK masses, fh RL is the EWSB mass matrix given in Eq. (55), and m p h ys 
is the diagonal matrix of physical masses (chosen to be real and positive). 

Also, from Eq. (56) one can obtain h—ri—ijj^—ip^ contact interactions involving the Higgs 



field, h. In the mass eigenbasis, these are controlled by the dimensionless couplings 



y-RL 
^ijk 



L Irak 



(57) 



with X$ 



^ikj ■ 



For self-annihilation processes, there are also relevant terms involving two radion modes 
and a fermion pair. Following Ref. |18j . we write these as 



00 00 

hi 2 =0jk=0 r 



14 



Some representative values were given in Tables 1 and 2 of that reference. 
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and 

oo oo 

= -£ E^K-^l^ + H . (58) 

*l*2=0jfc=0 r 

For a Higgs field localized on the IR branes, one has X?f- k = 4(Y" 5D /L) [l + (-l)»'i+*2+i+ fe ] x 
Fi(L)Fj(L)f( R (L)f$ L (L), where Fi is the wavefunction. In the same way as for the single 
radion interactions, the coefficients vX^j k , after rotation to the mass eigenbasis, can be 
assumed to be part of the m^jk a bove. 

The Feynman rules arising from Eq. (54) were given in Ref. [18J. We find it convenient 
to express various decay rates and cross sections involving the radion/KK-radion in terms 
of the dimensionless quantities 

M RL ■ ., M LR - ., 

G^., njh - n\ x , G u ihJk = n! x , (59 ) 

where rrij is the physical mass of the j-th KK fermion excitation, and we defined 

M ™4njk = g ('•"•'/m':,. ,/ + Cng^jkmk) ~ 'I n "I , jh , (60) 

with c n = (-3) n (2n/3 - 1) /n! and rf n = (-4) n -7n!. One also has M^\ jk = .\//;'.';, 
In the limit that EWSB is neglected, one has 



rn 



il—injk 



H L ..i nj k mj + 2g* R ln jk m k + rrv^jk , (61) 



where Ta^f n ^ k is a contribution coming from the Yukawa coupling involving the stabilizing 
scalar field, which is in general numerically negligible (the expression is given in [T8]). In 
this case, the g^ L in j k and g RR in j k are simply overlap integrals of n radion modes and two 
fermion wavefunctions. 



B Formulas for processes involving the Radion Tower 

In the following, we use extensively the notation 



rrii , , 

ejj = — . 62 
rrij 



Also, in the following formulas, the Gf- k , G^ jk , G R k and GfM are as defined in Eqs. (59) 



ijki ijki iijk^~ iijk 

ijk ^""^ i/ijA 



and (|60|), and the gf jk and g RR are as defined in Eq. (|54|). 



B.l Self- Annihilations 

For annihilation into fermion pairs, we are only interested in the process rV — > ff, but we 
give the result for the slightly more general case r^i — > f k f k . As shown in Fig. [TJ there are 
two types of diagrams: those involving a "heavy" / J exchange in the t- and ^-channels, and 
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a contact interaction. The contribution due to the t- and w-channel exchange can be written 

as 

N r mf 



r t+u 



2 \3/2 



~ \^ J " c ""i f^eff \ z (i _ 2 V 
W r,r^/'/' ~ l 6vr A 4 " 



(63) 



where we defined 



/^YCfT 

^* ijk 



^ 2 , I Ar <RL r LR 3 

t j,i t k,i > ^ i ^ J ijk (J ijk t j,i 



(l + e j,i _ e fc,i) ( 9 (fi^fc) + feffe) e *.» + [ 2 #i^A^ijf - 9ijk9ijk + ^9ijk^fjk\ e i,i 



,RR\ 2 



(64) 



and iV c = 3 (iV c = 1) for quarks (leptons), while the e's are the mass ratios defined in Eq. ( |62| , 
i.e. we normalize with respect to the mass. Eq. (63) holds in the deep non-relativistic 
limit (v is the relative velocity of the initial state). 



The non-relativistic contribution due to the contact interaction is 

, 2 



j. Contact 
riri^rf k f k 



N c m 2 k 



16tt Af. 



c2 

~k,i 



while the crossed-term is given by 



vcr 



Crossed 



EN C mjm k ff , RL nLR \ ( , 2 \3/2 
^ 4 u ijk \yiikk "T u iikk) \ l e k,i) 



(66) 



Annihilation into the KK-parity even Higgs degrees of freedom (see Fig. |2j), gives (see [18] 
for the relevant Feynman rules) 



N H Fi{Lf m 2 



1 - e 2 

1 t h,i 



4vr A$ (1 + e 2 - e 2 ) 2 



[e^(2 + 3e^)-ll e 2.-21e 4 M ] 2 , (67) 



where en,i = rnn/rrii and 6/^ = rrih/mi, with mu the mass of the "inert" Higgs doublet and 
mh the mass of the SM-like Higgs. Nh = 4 takes into account the 4 real Higgs d.o.f. in 
the limit that is large. The Fi(L) are the KK-radion wavefunctions evaluated on the 
IR boundaries, which are order one numbers, and in particular, F r /(L) « 1 for the first 
KK-radion mode. 

We also quote the tree-level annihilation cross section into identical massless gauge 
bosons, via a massive KK gauge boson exchange (of mass mj) in the t- and u-channels: 



(68) 



where £jj = mi/mj, and g^o is the overlap integral of the wavefunctions for a radion, a 
massive gauge KK mode and a massless gauge boson. This dimensionless coupling is of 
order 1 / y/2k eS L. 



37 



B.2 Co- Annihilations 

For coannihilations with the NLKP, involving gluons in the final state, we are interested in 
the process r't' — > gt, but we give the result for the more general case r t P — > gt k . Including 
the s- and t-channel diagrams, as well as a quartic contact interaction that is required by 
gauge invariance (see Ref. [IE]), we have: 



{N 2 C - IK 
16iV c A2 



1 + e 



J:' 



Cfc, 



(69) 



which holds in the ultra non-relativistic limit (v is the relative velocity of the initial state). 



The e's are the mass ratios defined in Eq. (62), i.e. we normalize with respect to the rj mass. 



There are also channels with a Higgs in the final state. We are interested in r't' — > ht, 
but we quote the more general case (in the ultra non-relativistic limit, and neglecting a 
subdominant diagram with a s-channel t k ): 



1 



ver, 



riV^rht k 



1 + e jti + e h . 



:i+e 



3,1 



'k,i. 



647rA? 



X 



[(*$) 2 + (*$) 2 ] [(l + ^) 2 + 4,-<J +^X^e kji (l + e j , i ) \ ,(70) 



where X-j% and X^ k were defined in Eq. (57). The e's are the mass ratios defined in Eq. (62), 
and we also defined = nth/mi, where rrih is the Higgs mass (i.e. we normalize all the 
masses to the r,- mass). 



B.3 p Decays 

The production of r» in /■?' decays is governed by 



1 P^nf k 



rrii 



16vrA2 
1 

x 



[1 - (e k ,j + Qj) 2 ] t 1 - - ^ 



J ( ' 



(71) 



where the e's are the mass ratios defined in Eq. (62), i.e. we normalize with respect to the 
mass of the decaying particle. 



B.4 Vj Decays 

The production of r, in V-[ decays, when the final gauge boson has mass m k ^ 0, is governed 
by 



-pCM 



1927im 2 k A 2 r 



1 - (e ktj + e itj y] [1 - (e feJ - e^-) 2 ] (l + e^. - e 2 ,) 2 + 8e 2 fej . ,(72) 
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where 

-<eff,V _ 



G$ V = 4^ fcCfcj - + g ijk (1 + e\ 3 - e^) , (73) 
and 



Qijk 



= ^ jjy Fifyfv , <S = ^ ' 2 , (74) 



with fl{y) the wavefunction corresponding to V<? (which is related to the V-> wavefunction 
by fl = dyfy/rrij). 

When nik = 0, we have instead 

2 3 

~ 167rA 2 ^ C mJ > 



where ^-fe is suppressed by order 1 / ^/2k e ^L due to the massless wavefunction. 
B.5 KK-parity odd Higgs Decays 

The states in the KK-parity odd Higgs doublet can decay into r' plus a KK-parity even Higgs. 
Parameterizing the KK-parity odd and even Higgs doublets as H = {H + , -^(h- +ia)} and 

H + = {G + ,v + ^(h + + iG )}, respectively, where v = 174 GeV, h+ is the SM-like Higgs 
and G°, G ± are the would-be Nambu-Goldstone bosons, we have 

F CM _ F r'( L ) 2m L h - 2 2 n2 



x y/ [l - (e h+ , h _ + e r , :h _f] [I - (e h+ , h _ - e r ,, h _f] (76) 

for the CP-even Higgs, and 

f cm _ F r'{L) 2 m 3 a , 2 , 3 f cm _ F r/ (L) 2 m 3 H± , 2 \3,__. 

for the CP-odd and charged Higgses. Here F r i(L) rs 1. 
B.6 Production of through Scattering 

The rj production cross section (times flux) via the process gqi — > (which includes s- 
and t-channel diagrams, as well as contact interactions required by gauge invariance) can be 
written as 

_a^_ m| 6(z,Ae,e ktj ) 

{[(0 2 + (OVi(^ Ae >^) + G S G Sf^(^A6,6 fej )} , (78) 
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where 



0(z,Ae,e Kj ) 
Fi(z,Ae,e k ,j) 



(1 + Aezy-AeljZ 



(79) 



tan 



1+Aez 



4 \1 + 2Aez+ (l + 2Ae + 2Ae 2 ) z 2 l n x 

9{z,Ae,e k)j ) 

- [3 + (2 + 7Ae)z + (3 + 2Ae)z 2 - Ae z 3 ] , 



10) 



F 2 (z,Ae,e K 



32e k j z 



[1 + Aez)- 



tan 1 



9(z,Ae,e M ) 
1+Aez 



9(z,Ae,e k% 



Here z = m|/s, with y/s the CM energy, Ae = — efj, and the e's are the mass ratios 
defined in Eq. (62), i.e. we normalize with respect to the incident P mass. The flux factor is 



given in the CM frame by F = y/s\p\, where \p\ is the momentum of any of the initial state 
particles, and replaces the relative velocity appropriate for the non-relativistic scattering of 
massive particles. 



The singularity at z — 1 in Eq. ( 78 ) arises from the s-channel diagram, and is regularized 



by the q 3 width, IY For instance, in the non-relativistic regime, z ~ 1, we have 



Fa 



a* m n 



gq^riq k 



8N C ( s 



m 2 ) 2 + m 2 r 2 



tan 



(l,Ae,e fc|j ) 
1+Ae 



e(l,Ae,e k j) 



x0(l,Ae,e feJ ) \[(G^ k ) 2 + (G^) 2 ](l + Ae)+AG5 L k G^e k A . (82 



In the ultra-relativistic regime, Eq. (|78) reduces to 



Fa. 



2 log 



mi 



(83) 



B.7 Radion Decays 



We also summarize the partial decay widths for the radion (KK-parity even) decays into 
fermion pairs, massive and massless gauge bosons, and the CP-even Higgs. These are used 



in Section |6.2| For the decays into fermion pairs, we have 

(^ooo) + (Gqm) (l - 2e 2 r ) - AGqqqGqqq e 2 r 1 , (84) 



T-*ff 



m r m) / 2 



where e/ jr = mf/m r , and a color factor N c should be included for quarks. 
For the decays into (indistinguishable) massive gauge bosons, we have 



F (L) 2 m 3 
32ttA2 



pmassive ^ 
1 r-i-VV ~ 



1-46^(1-46^ + 126^) 



15) 
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where ey, r = mv/m r , and Fq(L) ~ 1 is the wavefunction of the radion evaluated on the IR 
brane. Here we have assumed that the Higgs (hence EWSB) is IR-localized, and neglected a 
subdominant contribution from the bulk, present even for massless gauge bosons. The decay 
rate into massless gauge bosons is 

2 3 

■p massless 

T ^vv ~ ig^j, (86) 

where g oo ~ 1 / (2/c ff 2>) is the overlap integral for a radion and two massless gauge boson 
wavefunctions (the latter being flat). 

For the decays into the CP-even Higgs, h + , we have 

where = m^/mr, and a factor of 1/2 for identical particles has been included. A similar 
expression applies for the decay into the KK-parity odd Higgses, provided these channels are 
kinematically open. 
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